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A Survey of I. N. Vekua’s Theory of Elliptic Partial 
Differential Equations With Analytic Coefficients 
By PETER HENRICI, Los Angeles, California, USA?) 


In a series of papers which have appeared since 1937 and which were sum- 
marized in a book [14]?), the Russian mathematician I. N. VeKvA has developed 
a theory of linear elliptic partial differential equations with analytic coef- 
ficients and with two independent variables in a way that generalizes some 
aspects of the classical theory of holomorphic functions. VEKuA’s theory is of 
interest to the applied mathematician because, among other things, it provides 
him with an algorithm for the construction of so-called complete systems of 
solutions which in many-cases can be carried through in closed analytic form. 
In spite of their interest, VEKUA’s developments have been noticed but little 
by western mathematicians, and his book has not been translated. 

The present article is intended as an introduction to VeKua’s theory. 
We shall deal here mainly with his results concerning the analytic continuation 
of the solutions and their representation by holomorphic functions. VEKUA’s 
theory of the boundary value problems, which is intimately connected with 
these results, will be presented in a later report. 

In order to stress the practical usefulness of the theory, we have included 
in this paper a number of examples of equations which can be treated explicitly. 
Certain of these examples were not given by VEKUA. We also took the 
opportunity of including some other original material, such as the sections on 
singular equations (§ 4.3) and addition theorems (§ 4.4), and the entire para- 
graph on regular and singular elliptic Cauchy problems (§ 5). 

This article originated in a series of lectures delivered by the author at the 
Swiss Federal Institute of Technology during the spring term of 1956. The 
author is indebted to Professor E. STIEFEL, Director of the Institute for 
Applied Mathematics, for his advice and encouragement during the preparation 
of this report. . 

Notations 


Gothic capitals (6, §, D, ...) denote domains of the complex plane or of 


the real (x, y)-plane. 
0® denotes the boundary of ©. 


1) University of California. 
2) Numbers in brackets refer to References, page 202. 
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In the space K” of n complex variables z; (6. 1) 2, ~255 )y A Oa ee G,,) 
denotes the ‘cylinder domain’ obtained by letting each variable z, vary inde- 
pendently from the others in &,. 

For easier distinction between analytic functions of one and of several 
variables, regular analytic functions of one complex variable are called holo- 
morphic. 

Partial differentiations are mostly denoted by subscripts: 

odes iE 2 eas 2 ete 

A sequence of functions or an infinite series will be called uniformly convergent 

in a domain &, if the convergence is uniform in every closed subdomain of ©. 


1. Introduction 


In this section we shall lay down a few basic definitions and state the main 
theorem of VEKUA’s theory. 


1.1. Three Classes of Functions 


Let © be a domain of the (x, y)-plane. We consider three classes of complex- 
valued functions defined in G. We shall denote these classes by I(G), I1(G) 
and III (G) or, if there is no reason for ambiguity, by I, II and III. 

The class I(©) consists of all functions which are twice continuously differ- 
entiable in ©. Any functions of this class will be called regular in ©. 

The class II (©) consists of those functions f(x, y) which depend analytically 
on the two real variables x and y in ©. Functions of this class will be called 
real-analytic in ©. 

In order to define the class III (G), we recall some basic facts of the theory 
of functions of several complex variables. According to this theory there cor- 
responds to every /(x, y) in II(G) a function, called the analytic continuation 
of /(x, y), which is complex-analytic in a domain D of the space K? of the two 
complex variables x and y, containing ©, and which coincides with f(x, y) for 
(x, vy) € ©. We shall denote the analytic continuation of f(x, y) again by the 
symbol f(x, y). In general, there is nothing known about the size of D; the 
analytic continuation exists only ‘in the small’. 

Let us now introduce in K? two new variables z and z* by the relations 


Z=H+ty, st=x— ty. (ly 


The variables z and z* are conjugate complex if and only if x and y are real. 
Let f(x, y) € II (G); we put 


F(z, 2*) = f(x, y) . (2) 
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According to the above, every point z)= x9 + 7 479 € © possesses a neighborhood 
© (z9) such that F(z, 2*) is analytic in the cylinder domain [6 (z,), ©(z,)]. 
The class III (6), then, consists of all functions of class I1(G) for which 


G(z,) 0G forall 2,€¢6. 
This may also be stated thus: 


Definition 1.1. A function f(x, y) is said to belong to class II1(G), if the 
analytic continuation of the function 


exists and ts regular in the cylinder domain (@, ©). 

Thus, in a sense, for the functions of class III (©) the analytic continuation 
exists ‘in the large’. 

Example. Let © be simply connected, and let u(x, y) be a real harmonic 
function in ©. We assert that u(x, y) € III (G). 

Proof. It is known that there exists a function /(z) of one complex variable, 
holomorphic in ©, such that 


u(x, y) = Re f(z) . (1.3) 
We define 
f(z) = F(z) . : (ea) 


Clearly, f(z) is holomorphic for z € 6. Let us now consider the function 


Ute, 2*) = > [H@) + fe). 


U(z, z*) is analytic for (z, z*) € (G, G). By (1.3) and (1.4), 


Ue, 2) = > [#le) + 1] = u(x, 9) - 


Hence U(z, z*) represents the analytic continuation of (x, y) into the cylinder 


domain (6, G). 


1.2. The Differential Equation Considered 


We consider the linear elliptic partial differential equation of the second 
order and with two independent variables for an unknown function u(x, y) in 
its normal form 


e(u) = Au+au,+ bu,+cu=0 (e) 
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(Au = U4, + U9). The functions a= a(x, y), b= d(x, y) and c=c(x, y) are 
called the coefficients of the differential equation. The conditions to be satis- 
fied by them will be stated below. 

Besides equation (e) we shall also consider its adjoint equation, defined by 


é(u) = Au —(au),—(bu),+ou=Au+dumtbu,tiu=0. (8) 
It is easy to verify that é(u) = e(u). 
VEKUA’s theory can also be extended to elliptic equations of higher order 
and to systems of elliptic equations, but we shall not discuss these generaliza- 
tions here. 


1.3. The Theorem of Vekua 
The following is a classical result due to PICARD: 


Theorem 1.1. For arbitrary ©, if the coefficients of (e) are in II (G), then 
every solution in I (G) ts also in I1(G). 
The theorem of VEKUA can be stated analogously thus: 


Theorem 1.2 (main theorem). If © 1s simply connected and 1f the coefficients 
of (e) ave in III (G), then every solution in I (G) ts also in III (G). 


Table 1 
The Theorems of Picard and Vekua 


| PICARD | VEKUA 
If domain arbitrary simply connected 
coefficients TT UE 
solution if I 
then solution | II | III 


We can formulate VEKUA’s theorem more smoothly if we introduce the 
concept of a fundamental domain. 


Definition 1.2. A domain © 1s called a fundamental domain of equation (e), 
if it 1s simply connected and if the coefficients of (e) belong to III (6). 

We note that if the coefficients of the equation are polynomials or entire 
functions, then the whole plane is a fundamental domain. For any equation 
with analytic coefficients, every sufficiently small domain is fundamental. Also, 
any simply connected subdomain of a fundamental domain is itself funda- 
mental. A fundamental domain of (e) is also a fundamental domain of (é). 

We now can state VEKUA’s theorem thus: 
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Theorem 1.2. If © is a fundamental domain of (e), then every regular 
solution belongs to III (G). 

Evidently this theorem is a generalization of the statement of the example 
of section 1.1. The same example shows that the theorem can no longer remain 
true if © is multiply connected. In fact, the function 


u(x, y) = log (x? + y?) 


is harmonic in every annular domain surrounding the origin. But the 
analytic continuation of u(x, y), 


U(z, 2*) = logz + logz* , 
is no longer analytic in (, ©). 

For simply connected domains, the statements of both Theorem 1.1 and 
Theorem 1.2 can be characterized by saying that ‘the behavior of the solutions 
imitates the behavior of the coefficients’. In the Vekua case in particular, any 
representation of the solutions in ITI (G) will also hold for all regular solutions, 
provided that the coefficients are in III(G). Such general representations of 
the solutions will be derived in sections 4 and 5. The proof of the main theorem 
will be given in section 3. In section 2 we shall introduce an important auxiliary 
function. 


2. The Riemann Function 


In this section we shall introduce a special solution of (e) which will be of 
use in all later developments. This solution depends only on the differential 
equation and not on the particular domain considered. It is called the Riemann 
function of e(w) = 0 because of its formal analogy with the functions thus 
called in RIEMANN’S classical method of integrating hyperbolic differential 
equations. It is defined as solution of an integral equation of the Volterra type 
and can always be constructed by the method of successive approximations; 
for some equations it can be expressed explicitly in terms of known functions. 


2.1. Volterra Integral Equations in the Complex Domain 


We first consider the case of an equation with an unknown function depend- 
ing on one variable. 

Let © be a simply connected domain, let f(z) be holomorphic in ©, and let 
K(z, t) be holomorphic in the cylinder domain (6, ©). Let 29 be a fixed point 
in ©. We consider the functional equation 


w(e) — | K(@,#) w(t) dt = fle) (2.1) 
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and call any function w(z) which is holomorphic in © and satisfies (2.1) a solu- 
tion of the equation. 
Theorem 2.1. The integral equation (2.1) possesses exactly one solution. 
Proof. For }(z) = 0 the Taylor expansion of w(z) at z = 2 vanishes identi- 
cally. This shows the uniqueness of the solution. To prove the existence of a 
solution, we construct a sequence of functions w”)(z) (n = 0, 1, ...) as follows: 


w)(z) = f(z); w(z) = f(z) +/Ke t) wl OG) dt = 1, 2, -) ee 


It is evident that if 
lim w")(z) = w(z) (2.3) 


exists uniformly in ©, then w(z) is a solution of (2.1). To demonstrate the 
uniform convergence, we observe that 


w(2) = fz) + is I (z, t) f(t) dt, (2.4) 
where . 
P(g, 2) = KM(z, t) + K(z, t) 4. + KO, 2) (2.5) 
and 
Kz, t) = K(z,t), K@*%(z,2) = a K(z, s) Ks, t) ds. (2.6) 
t 


The sum on the right of (2.5) is easily seen to converge uniformly in (©, G) by 
the principle of dominated convergence, and this establishes (2.3). 
Defining 


1A? y KOZ ae (2.7) 


we can write the solution of (2.1) in the form 
w(2) =f) + (Pe) 10 ae. (2.8) 


The function J'(z, ), which depends neither on f(z) nor on Z, but only on the 
kernel K(z, ¢) of the integral equation, will be called the resolving kernel of the 
equation. 

We add three corollaries to the above results. 


Corollary 2.1. The solution w(z) of (2.1) is a holomorphic function of z for 
2 € G. 
This is an immediate consequence of (2.8). 
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Corollary 2.2. If for every (z, t) € (G, ©) the kernel is a holomorphic function 
of a parameter € in a fixed domain §, then the resolving kernel I’ = I(z, t, 6) 
1s analytic for (z, t, C) € (6,G, §). 

The uniform convergence of J"(z, t) = I"(z, t, ¢) can be established as above. — 
If also the right hand side of (2.1) is a holomorphic function of ¢ in §, then the 
solution w = w(z, ¢) is analytic in (6, §). 


Corollary 2.3. For (z, t) € (©, G), the identity 


T(z, t) — / nis s) I(s, t) ds = K(z, t) (2.9) 


holds. 
From (2.5) and (2.6) we have 


(zt) — K(z, s) I'"-9(s, t) ds = K(z,t) . 
t 


The required result follows by letting 7 > oo. 

The main task of this section consists in establishing analogous results for 
the following more complicated equation, which involves integrations with 
respect to two variables: 


w(z, 2*) [Kies *, t) w(t, 2*) dt — - [i a ear(een*) ai 
(2.10) 


a* 


fi dt vi Ke: 2* t, t*) w(t, t*) dt = F(z, z*). 


Here the functions K,, K,, K, and F are assumed to be analytic for 


(2, 2*; t, t*) € (6, H; GH); 


© and § are two simply connected domains. The points z) and z* are fixed in 
G and §, respectively. 

We shall prove for equation (2.10) the following analogues of Theorem 2.1 
and of Corollary 2.1: 


Theorem 2.2. Equation (2.10) possesses a unique solution w(z, z*) which ts 
analytic in (G, §). 


Corollary 2.4. The solution w(z, z*) depends analytically on (Zo, 2) in the 
domain (G, §). 
We shall not need the analogues of the Corollaries 2.2 and 2.3. 
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Proof. The methods which were used to deal with the simple equation (2.1) 
can be carried over without change to the case K, = K, = 0 of equation (2.10). 
It is thus easy to show that the ‘reduced’ equation 


wy sf dt li K,(z, 2%, t, t*) wold, t*) di® = F(z, 2*) , (2.11) 


where K, satisfies the same conditions as K in (2.10), possesses a unique 
solution w,(z, z*) which is analytic in (G, §). Also Corollary 2.4 is easily seen to 
hold in this case. 

In order to deal with the complete equation (2.10), we consider the functions 
K,(z, 2*, t) and K,(z*, z, #*) as kernels of simple integral equations of the type 
(2.1) with the respective domains © and §. The kernel K, (K,) is an analytic 
function of the parameter z* (z) in the domain § (G). According to the Corol- 
laries 2.2 and 2.3, the resolving kernels I’, = I,(z, z*, ¢) and I, = I’,(z*, z, ¢*) 
are analytic functions in (G, §, G) and (H, G, H), respectively, and satisfy 
the relations 


0 — [Biles LAS 2850) GS Ial eee ee 
(2.12) 
T(z, 2, 08) — / K,(z*, z, s*) D'g(s*, z, 8) ds* = Ko(2*, 2, t*) . 


t* 


We now seek the solutions w/(z, 2*) of the complete equation (2.10) in the 
form 


(2, 2") = (2, 2*) + {Tie 2*, t) Wot, 27) at mee (2%; 2,0*) wele et eat™ : 


(2.13) 


where w,(z, 2*)-is a new unknown function. Introducing this into (2.10) we 
find after some cancelling, using (2.12), a reduced equation for w,(z, z*) of the 
form (2.11), where 


I Cc dale oi hd Wee SLR ee ee 
+ K,(2, 2, 1) P(2*, 2, ) + K,(2*, 2, 2) Ty 2%, 0) 


+ [Ke K(z, 2*, s, t*) I'y(s, ¢*, 4) ds (2.14) 


+ [Ke at, tS‘) se*) case 


te 
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Since this equation possesses a unique solution Wy(z, 2*), which is analytic in 
(G, ), we can reverse our argument and conclude that w(z, 2*) as defined by 
(2.13) is a function — and also the only function—which solves (2.10). This proves 
Theorem 2.2. Furthermore, since ws(z,z*) depends analytically also on (zo, 2%) 


in (©, §), we can again by the representation (2.13) verify the assertion of 
Corollary 2.4. 


2.2. Definition and Fundamental Properties of the Riemann Function 


Upon substituting the variables 
Z=xKKty, = x4—-—7y (1:1) 
in (e), one obtains for the function 


gt — * 
vier SE 


2 
4 


the differential equation 


SAA BAD, HEU S CULO, (E) 
where 
A(z, 2*) y=ifa (x, y) +7 0(x, y)}, 
i 
Biz, 2*) == {4 a(x, y) —10(x, y)}, | (2. 5) 
1 
C2") == ¢ C(x y) 


Similarly the adjoint equation (é) transforms into 
E(U) = U,., — (A U),—(BU),+ CU= Pips AW 8 Uerk Cul Oeb(k) 


Thus, the two operations of passing to the adjoint equation and to the variables 


(1.1) can be interchanged. . 
If & is a fundamental domain of (e), the functions A, B, and C are analytic 


in the cylinder domain (6, 6). Therefore for (z, 2*, t, ¢*) € (G, 6, 6, 6) the 
integral equation 


V (Gz, 2") ines 2z*) V(s, 2*) ds vA ot} aes Sthas* ae 


+ [ufo (s, s*) V(s, s*) ds*¥ =1 


is of the type studied in section 2.1. According to Theorem 2.2 and Corollary — 
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2.4, it possesses a unique solution which is an analytic function of (z, 2*, Z, ¢*) 
in the quadruple cylinder domain (©, ©, ©, G). 


Definition 2.1. The solution of (2.16) is called the Riemann function of 
equation (e) and will be denoted by R(z, 2*, t, t*). 

We shall state some fundamental properties of the Riemann function. 

Theorem 2.3. With respect to its first two arguments the Riemann function. 
of (e) satisfies the adjoint equation (é). 

Application of the operator 0?/0z 0z* to the integral equation (2.16) yields 
(E), the complex form of (é). 


Theorem 2.4. For (z, 2*, t, t*) € (G, 6, 6, 6), the following relations hold: 


WAG bey sae oA a bs : (2.17) 
Rylent* it) Sav OG, Peace pa (2.18a) 
Ro(G,-2*; 60") =e Aha") ie (2.18b) 
Ralape" yt, S70) = — BD Ge eee oe eee (2.18c) 
Ky 2) 2% 12, 0) = A (2 ee ee (2.18d) 


Relation (2.17) is evident. To prove (2.18a) and (2.18c), we put in (2.16) 
2° ==)" and get 


R(z, *, t, 1%) — | Bos t*) R(s, t*, t, ) ds =1. (2.19) 


The required relations follow by differentiation with respect to z and ¢, respec- 
tively. The relations (2.18b) and (2.18d) are proved similarly. 

It is easy to see that the properties expressed by Theorem 2.3 and by the 
relations (2.17), (2.18a) and (2.18b) characterize the Riemann function com- 
pletely. ; 

The following theorem concerns the Riemann function - the adjoint 
equation. 


Theorem 2.5. J} Riz, z*,t,t*) 1s the Riemann function of the adjoint 
equation (é), then _ 
Rl, 25 0 te ee ee ees (2.20) 


We can thus obtain the Riemann function of the adjoint equation simply 
by interchanging the first and the second pair of arguments. 
Proof. Since the relations (2.18a) and (2.18b) for R amount to (2.18c) and 


(2.18d) for R, we have only to show that R(é, ¢*, z, z*) as a function of z and 
2* satisfies E(R) = 0. 
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Let U = U(s, s*) be an arbitrary function analytic in (6, 6), and let 
R= R(s, s*, z, 2*). Then the following identityis true 


(UR) — RE(U) =[U(R,— A R)], + [U(R, — B R)]q. (2.21) 


Integrating with respect to s and s* between the limits ¢, z, and ¢*, z*, respec- 
tively, where (t, ¢*) is a fixed point in (6, G), we obtain 


Pia at) = UG, t) RE, es, 2?) 


ip | R(s, t*, 2, 2*) [U,(s, t*) + B(s, #*) U(s, t*)] ds 


. (2. 22) 
+ / Rit, s*, z, 2*) (U(t, s*) + Ale, s*) Ul, s*)] ds* 
a 
+ | ds | R(s, s*, 2, 2*) E(U(s, s*)) ds*. 


t t* 


Let now U(z, z*) = R(t, t*, z, z*). In virtue of (2.17), (2.18c) and (2.18d) we 
then obtain é < 


[as [ R(s, s*, 2, 28) E(RU, t%, s, s*)) ds* = 0. (2.23) 
. eh 


From this the assertion follows, since all functions involved are analytic. 


2.3. Transformation of Variables and the Riemann Function 


The Riemann function is known explicitly only for very few differential 
equations (see Table 2). Since other equations can sometimes be reduced to one 
of these standard types by transformations of the dependent or independent 
variables, it is important to know how the Riemann function behaves under 
such transformations. We shall state here two theorems concerning this be- 
havior. The proofs are purely formal and are therefore omitted. 

Let e(w) = 0 have the fundamental domain 6, and let E(U) = 0 be the 
complex form of e(w) = 0. By introducing the new dependent variable 


ae z*) Fal ag an 


U(z, g*) 
where ®(z, z*) is analytic in (©, ©) we obtain a new differential equation 
E(0) =U, +AU, +BU,4+CU=0, (2.25) 
where 
E(@). (2.26) 


ol 


A=A+(log®),, B= B+ (log%),, C= 


ZAMP 
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Theorem 2.6. If O + 0 for (z, z*) € (G, 6), the equation E(U) = 0 possesses 
again the fundamental domain ©, and its Riemann function is given by 
O(z, 2*) 


Bid A ale A ta Oe b(t, ) era a: Edt) (2222) 


Let ©’ be a simply connected domain which is mapped conformally onto © 
by the transformation é ; 
afc}, = ¢ € B", 


and let 
ae = f(f*), C*e 6’, 


where /(¢) = f(#). Let E(U) = 0 be the complex form of a differential equation 
e(u) = 0 with fundamental domain © as above. Then the function 


satisfies the equation 
BP i. AU + BU, CU" 0 (2.28) 


with the coefficients 


, = yi; \ d, 
BYE, o*) = BAO), Ae) AE, (2.29) 
‘ Tears eee 
CUE, 9) = C(H0), Hon) 
Theorem 2.7. The equation E'(U') = 0 has the fundamental domain 6’, 
and its Riemann function is given by ; 


BEG eS, ) = RU (2). HE*), Seer) (2. 30) 


3. The Fundamental Solution 


In this section we shall prove VEKUA’s Theorem 1.2. As in the usual proof of 
the classical Picard theorem, the proof is accomplished by expressing the 
solution in terms of the boundary values, using a so-called fundamental 
solution; the stronger result follows by making fuller use of the analyticity of 
the fundamental solution in the large. 


3.1, Construction of a Fundamental Solution 


Definition 3.1. Let © be a fundamental domain of e(u) = 0, and let Q = (&, ) 
be a point of ©. Any solution of e(u) = 0 which in © — Q behaves like 


e(x, y) logr + go(x, 9) (r=V (x — 8)? + (vy --0)?), 
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where the functions g(x, y) and g(x, y) are regular in © and where g(é, y) + 0, 
is called a fundamental solution of e(u) = 0. The point Q 1s called the pole of the 
fundamental solution. If g(E, 4) = 1/20, the fundamental solution is called nor- 
malized. 

By adding any regular solution to a fundamental solution or by multiplying 
it by a nonvanishing constant one obtains new fundamental solutions. The 
fundamental solution of a given equation and with a given pole is thus not 
uniquely determined. 

The actual construction of a fundamental solution is carried out with the 
help of the following lemma. 


Lemma 3.1. Let the closure of the simply connected domain © le in a fun- 
damental domain of e(u) = 0, let the functions y(z) and p*(z*) be holomorphic in 


G and G, respectively, and let (¢, €*) be in (G + 0G, G+ 06). Then, provided 
that the integrals exist, the two expressions 


U(z, 2*) = lz) Rlz, C*, 2, 2%) — ii os) Ry(s, C*, 2, 2*) ds , (3.1) 
€ 


ok 


Ue se" =O (27) (Gage ee) = gp (s*) Ke (€, 3%, 2, 27 gst Kaan 
c* 
are solutions of E(U) = 0. 
The proof is by straightforward verification, using Theorem 2.5 and the 
relations (2.18). 


Let us now assume that ¢ € 06, ¢* € 0G. Putting in (SL }, (32) 


plz) =log(z—¢), p*(z*) = log (2* — ¢*) , 


we obtain after simple calculations the two solutions 
al 
U2, g* hie HCC a2") log ig — C) — [logs ee ta + $(z— Oy C*; z, z*) Ses 
; ' 


OM, at) RC" 7s, ot) doe lee — Ce) 


0 
— flogs* aa# RC, CPs (Fee a) as*. 
We write 


A(z, 2, £, £*) = <o_{ Ulz, 2*) + U*(z, 2*)}. (3.3) 
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Evidently this function can be written in the form 
Soe 1 = - ‘ 
A(z, 2%. C; &*) =a 4 R(¢, G*, Zz, 2*) log | (z re a (2* - é*)| sae Ay(z, at, By Gr) , 


where the function A,(z, 2*, €¢*) is analytic in (6, 6, 6, G). The following 
statement is now obvious: 


Theorem 3.1. Let © be a fundamental domain of (e) and let § be a simply 
connected domain such that § + 0$ C G. Then for z€ d§, 2* € ds, the function 
A(z, 2*, &, &*) and its derivatives with respect to z and z* depend analytically on 
(¢, ¢*) in (, 8). 

For z* = z, if z moves around the point € along an arbitrary closed contour, 
the expression 


b| Fed 


log (z — ¢) @ — ¢) = logr 


returns to its initial value. Also, G(¢, Bee ¢) = 1 in virtue of (2.17). We there- 
fore have established 


Theorem 3.2. Let Q = (€, 7) © G. The function 
A(x, ¥p 6; 4) a Ale +4 yy, x —4 ve, € +4, § — 1) (3.4) 


is a normalized fundamental solution, regular in © — Q, of e(u) = 0. 


3.2. Proof of Vekua’s Theorem 
The following generalization of GREEN’s formula is well known: 


Theorem 3.3. Let § be a (not necessarily simply connected) domain with 
piecewise smooth boundary 0, and let v be the interior normal of 0%. Let the 
coefficients of e(u) be in class 1(§), and let u.and v be two arbitrary functions 
in 1(§). Then the following identity holds: 


(6.5) 
— {w we —— = 40 [a cogty, x) 6cos(y, y)]h ds . 


~ Let now & be a fundamental domain of e(w) and let u(x, y) be an arbitrary 
solution belonging to I (G). In order to show that the function U(z, 2*) = u(x, y) 


‘can be continued analytically into (6, 6), it is sufficient to show that the con- 


tinuation is possible into every cylinder domain (D, D) where D+ OD CG. 
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Let D be a domain of the required kind, and let Q = (&,7) € D. According 
to Theorem 3.2, the adjoint differential equation é(w) = 0 possesses a normal- 


ized fundamental solution Ax, y, €,y) regular in © — Q and having its pole at 
Q. We now choose in GREEN’S generalized formula (3.5) § = D — &,, where 
R, denotes the disc |z— €| < @ (¢ =&+7y), and put 


v(x, y) = A(x, y, & 7) - 
Letting @ > 0 we obtain 


u(g,n) = [[u % — anol as, (3.6) 
0d 


where on the right w = u(x, y) and 


N(u) = “* + aucos(r, x) + bucos(y, y). (3-7) 


The integration is with respect to the point (x, y). 
Formula (3.6) expresses the value of w at an arbitrary interior point of D 
in terms of the values of w and the first derivatives of w on the boundary of D. 


We now substitute in (3.6) for Ax, y, €,n) the function A(z, Gul ay: The 
normal derivative dA/dy can be expressed as a linear combination of A, and Ay. 
We thus obtain for U(¢, ¢) = u(&, 7) an expression of the form 


U(6, 0) = | {A+ tert fede} del, (3.8) 
oe) 


where the functions hy f., and E depend only on the boundary point z. According 
to Theorem 3.1, A, A and Ay as functions of the last two arguments can be 


continued analytically into the cylinder domain (D, D). By integration with 
respect to the parameter z this property is preserved, and thus the same is true 


for U(C, é). This completes the proof of VEKUA’s theorem. 


3.3. Riemann’s Function and Fundamental Solution 


In section 3.1 we derived a fundamental solution from the Riemann function. 
In this section we shall show that it is also possible to obtain the Riemann 
function from an arbitrary fundamental solution. As a byproduct we shall find 
that the coefficient of the logarithm in the fundamental solution is uniquely 
determined. . 

Let G be a fundamental domain of e(u), let Q = (é, 7) € G and 


A(x, y, €,) = e(%, y, &, 9) logr + a(x, y, &,n) 
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be a fundamental solution of e(w) = 0, regular in © — Q, with the added 
condition that the functions @ and o belong to class III() with respect to 
(x, y) as well as (€, 7). 

We put 


P(z, 2%, 6 C*) = R(C, C*, 2, 2*) , (3.9) 


1. €., apart from the order of the two. pairs of variables the Riemann function is 
identical with the coefficient of the logarithm in a normalized fundamental solution. ., 

For the proof it suffices to show that R as function of (z, z*) fulfills E(R) = 0 
and that it satisfies the relations (2.17), (2.18c) and (2.18d). This will show 
that R is the Riemann function of the adjoint equation é(u) = 0; by Theorem 
2.5 this establishes the required result. 

Utilizing the fact that A(z, 2*, ¢, C*) = A(x, y, é,) satisfies the given 
differential equation, we find 


Bi) EEG SUES ae (ix ANE 


1 
2 ee C*) 


(3.10) 


2 {P,+ B P}+ regular function = 0. 
From the fact that an identically vanishing function can have no logarithmic 
singularities, we conclude that E(P) = 0. Since it also can have no pole-like sin- 
gularities, we obtain the relations 


Ps et ea Algo) Llane eG = 0, 
Pele oe, S) a Be 2) 2 eet 2) =H 05 
which are equivalent to (2.18c) and (2.18d). The proof of (2.17) is trivial. 
We might add that, contrary to Theorem 2.5, the function 2(&, 7, x, y) is 
in general not a fundamental solution of the adjoint equation. 
In concrete cases Theorem 3.4 can be useful for the purpose of representing 


the Riemann function explicitly. The reader is urged to construct the Riemann 
function of the equation of axisymmetric potentials, 


Au + 5 ty = 0; haeiesetl) 


by superposing the three-dimensional fundamental solution of the potential 
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equations, 
o(x, y, 9; p) =[(x — 62+ (y— cose)? + (y—gsing)"]" (0Sp<2zn), 


verifying that there is a logarithmic singularity at (x, y) = (&, 4), and extracting 
the coefficient of the logarithm. 


4. General Representations of the Solutions 


In this section we shall derive general representations of all solutions ot 
e(u) = 0 which are regular in a fundamental domain in terms of holomorphic 
functions. Similar representations have been given by VEKUA for solutions 
defined in multiply connected domains, and for solutions of equations of higher 
order and of systems of elliptic equations’). 


4.1. The Goursat Problem 


Let &© be a fundamental domain of e(w) = 0, and let (¢, €*) be a fixed point 


in (G, G). Let y(t) and g*(¢*) be two functions which are holomorphic in © 
and @, respectively, and satisfy the condition 


(o) = p*(c*) - (4.1) 


We shall study the following problem (which for more general equations, but 
only in the small, was first considered by GourRsarT): 

To find a solution U(z, z*) of E(U) = 0 which is analytic in (6, 6) and 
satisfies the conditions 


Ut, C*) =o), 1€G; Ul, =q*(t*), *CGC. (4.2) 


Theorem 4.1. The Goursat problem has a unique solution. 
Proof. The solution, if it exists, satisfies the equation 


i= apr ap wy = euet (4.3) 


where A, B, and C are the coefficients of E(U) = 0. Integrating (4.3) with 
respect to the first variable from ¢ to z and with respect to the second from 


3) Representations of solutions of elliptic equations in terms of holomorphic functions have 
also been given by S. BERGMAN [1]. From some points of view these representations are more 


complicated than those given here, and the conditions under which they are valid are not easily 
summarized. 
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¢* to z*, where (z, z*) is an arbitrary point of (G, 6), we find the identity 


z ok 


U(z, 2*) — | Be c*) Ut, f4).dt — [AC #) U(E, t*) dt* 
aig p (4-4) 
+ | dt | C(t, *) Ult, t*) dt* = D(z) + O*(z*) — OL) , 
where Sa 
D(z) = plz) — [ BY L*) pli) dt, O*(2*) = g*(e4) — [AC ptt) ae. (4.5) 


“ 
i 
S 


The integral equation (4.4) is of the type (2.10) and possesses by virtue of 
Theorem 2.2 a unique solution which is analytic in (©, 6). This shows that 
there is at most one solution of the Goursat problem. On the other hand, by 
defining U(z, z*) as the solution of (4.4), we can easily show that it satisfies the 
equation (4. 3) and the conditions (4.2). This proves Theorem 4.1. 

We now shall represent the solution in terms of definite integrals involving 
the Riemann function. If as before R(z, z*, ¢, #*) denotes the Riemann function 
of e(u) = O, it satisfies by Theorem 2.5 the integral equation 


R(s, s*, z, 2*) — [Be 2*) R(s, s*, t, 2*) dt — i A(z, t*) R(s, s*, 2, t*) dt* 
ath $ (4.6) 
+ [ ae/ Cl 1*) R(s, s*, t, f*) df= 1, 


which holds identically for (s, s*, z, z*) € (G, 6, G, G). Putting s* = ¢*, multi- 
plying by @’(s) and integrating with respect to s from ¢ to z, we find that the 
function 


U2, 2) = i @'(s) R(s, C*, 2, 2*) ds (4.7) 
¢ 


satisfies 
V(U™) = Bz) — B(¢) , (4.8) 


where V(U) denotes the operation on the left of (4.4). Similarly 
Uz, 24) = | OY (st) RUC, st, 2, 24) dot (4.9) 
c* 
can be shown to satisfy 
V (GS Ot 2") OF (6*).. (4.10) 


188 PETER HENRICI ZAMP 
Since obviously 3 Z . 
U®(z, 2*) = D(C) R(C, C*, 2, 2*) (4.11) 

fulfills 
V(U®) = @(2), (4.12) 


we see that the solution of (4.4) can be represented in the form 


CS eae U®) . 
We thus have proved 


Theorem 4.2. The solution of the Goursat problem can be represented in the 
form 


Ue, 2) DiC) RWG Cee) + [8 Rik cece ee 
1 (4.13) 


ee 7 D*'(t*) R(C, t*, 2, z*) dt*, 
(a0 
where D(z) and B*(z*) are given by (4.5). 
The importance of Theorem 4.2 rests on the fact that the representation 
(4.13) holds not only for the solutions of the Goursat problem but also for all 
solutions which are only supposed to be regular in the fundamental domain. 


Theorem 4.3. Let © be a fundamental domain of e(u) = 0 and let u(x, y) 
be a solution which 1s regular in ©. Then there exist two functions y(z) and o*(z*), 
holomorphic in © and ©, respectively, such that u(x, y) can be represented in the 
form (4.13) (with 2* =2z,z=x+7,y). 

Proof. By VeKua’s Theorem 1.2, u(x, y) can be continued analytically 
into (G, G), i.e., there exists a function U(z, z*) which is analytic in (6, 6) 
and satisfies ; : : 

U(ixtiy,x—ity)=ulx,y), (*%, y)E®: 
We define g(z) and *(z*) by 
P(e) = U(e.C* i Ea) = UC. (4.14) 


The functions g and g* are holomorphic in © and 6. Let us now solve the 
Goursat problem with the conditions (4.14). The solution can be represented 
in the form (4.13). In view of the uniqueness of the solution it must coincide 
with U(z, 2*) and, for z* = z, with u(x, y). 


4.2. Approximation of the Solutions 


Let © be a fundamental domain of e(u) and let QO = (¢, ¢*) bea fixed point 


in (6, G). Let the functions yz) and g*(z*) be holomorphic in & and 6, 
respectively. 
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Definition 4.1. The pair of functions 


P = [plz), v*(z*)] 
is called admissible (with respect to the point Q), if 


P(e) = p*(c*) . (4.15) 


It is evident that for fixed G and Q the set of alladmissible pairs forms a vector 
space. If p, = (@,, pf) and 2 = (2, p¥) are two admissible pairs and if « 
and f are two arbitrary complex numbers, we define the pair # = « fp, + B p; by 


p= (a pit B po, a Pf + B DF) . (4.16) 


If for a sequence of admissible pairs ~, = (~,, y*) (2 = 0, 1,...) there exists 
an admissible pair # = (p, p*) such that y, > p, p* > p* uniformly in G and 
(, respectively, we say that the sequence #, converges to ~, and we shall 


write 
Pn>p (n>oo). 


We may look at the expression on the right of (4.13) as an operator which 
maps the space of admissible pairs onto the space of regular solutions. We shall 
write accordingly 


u(x, y) = QP). (4.17) 


According to Theorem 4.3, the operator Q possesses an inverse. If u(x, y) is a 
regular solution, then 

2Q-Tu] = [U(s, c*), U(C, 2*)], (4.18) 
where 


piel Zz = 2k z— 2* ) 
Ute; 2) w ( Sat eeree iri ued Je 
Moreover, the following theorems are easy consequences of the represen- 
tation (4.13): 


Theorem 4.4. The operator Q is linear. 
This means that for any two admissible pairs f, and f, and for any two 
complex numbers «, 8, 


Qa py t+ B pe] = & Spi] + B QApe] - (4.19) 


Theorem 4.5. The operator Q is continuous in the following sense: If a 
sequence of admissible pairs p, (n = 0,1, ...) converges to an admissible pair p, 
then 

Qn) > 261 (4.20) 


uniformly in (G, ©). 
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As a consequence of these theorems, every method of approximating holo- 
morphic functions leads to a method for approximating solutions of e(u) = 0. 
For the sake of illustration we discuss the approximation by what we shall 
call complete systems of functions. 


Definition 4.2. Let © be a domain of the complex plane, and let the func- 
tions f(z) (n = 0,1,...) be holomorphic in G. We shall call the system { fn(z) } 
complete with respect to © if every function f(z) which is holomorphic in © admits 
a unique expansion of the form 


fle) = 3 ay fale) (4.21) 


which converges uniformly in ©. 

It is well known that there exist (infinitely many) complete systems of 
holomorphic functions with respect to every finitely connected domain ©. 
Some examples are listed in Table 3. 


Table 3 
Some Complete Systems of Holomorphic Functions 


Domain _ Functions f, (2) 
Invervorote moles (COM (ame iN aeons ze 
Interior of ellipse with fociatz =+ 1 . | Jacobipolynomials Pi%B) (z), in particular: 


Legendre polynomials P, (z). 


n 


Interior of strip Im'2)= ‘const ~~ 2. Hermite polynomials 4H, (z)+) 


Interior of parabola Re (— z)1/? = const. | Laguerre polynomials L%(z)}) 


1) These systems, however, are complete only with respect to a subclass of the class of 
functions which are holomorphic in ®. In order that f(z) be expandable in the form (4. 21), 


it must satisfy certain growth conditions (see [10] and [12]). 


Theorem 4.6. Let © be a fundamental domain of e(u). There exists a system 
of solutions Un(X%, y) (n= 0,1,...) of e(u) = 0 such that every solution in 1(6) 
admits a unique expansion, uniformly convergent in &, of the form 


u(x, ¥) = 3 ay ttyl, 9) (4.22) 

; n=0 
' Proof. Let Q = (¢, ¢*) be a fixed point in (6, 6), and let 7,,(2) (== 0, 1,2 
be a complete system of holomorphic functions with respect to G. The system 


f(z) = fn(Z) is complete with respect to ©. We define the pairs 


Po= (1,1), Pon =L[(z—C) fo-a(2), O1, 
Poni = (0, (@* — 0°) rate) = 1,2,...). ue 
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Clearly, the pairs #, are all admissible with respect to Q; moreover, the system 


of pairs ~, (1 = 0,1, ...) is complete in the sense that every admissible pair £ 
admits a unique expansion of the form 


P= Db) an Pa: (4.24) 
n=0 


We define 
u,(x, ¥) = Q[D,] - (4.25) 


Let now u(x, y) be an arbitrary regular solution. In virtue of Theorem 4.3 there 
exists an admissible pair p such that u = Q[p]. By (4.24) and (4.25), 


= al 5 an» P| = yy Ay Uy (x, y) , 
n=0 n=0 


which is an expansion of the desired type. The uniqueness of the expansion 
follows from the uniqueness of (4. 24). 


Definition 4.3. The system of solutions u,(x, y) (n= 0,1,...) of Theorem 
4.6 ts called a complete system of solutions of e(u) = 0 with respect to the domain ©. 

The complete system (4.25) is of course not fully determined by the differen- 
tial equation and the fundamental domain. It depends on the choice of the 
complete system of holomorphic functions and on the point (¢, ¢*). On the other 
hand, there exists one distinguished system of solutions which frequently can be 
constructed explicitly and which therefore is important from the point of view 
of applications. This is the system (4.25) where € = ¢* = 0 and 


Po ee (1,1) , Pon om (2", 0) ’ Pen-1 =i (0, Fes (4. 26) 


It can always be defined when the origin can be included in a fundamental 
domain. 


Definition 4.4. The system of solutions given by (4.26) (with ¢ = ¢* = 0) zs 
called the normal system of solutions of e(u) = 0. 
Some examples of standard systems are given in Table 4. 


Since for circles around the origin the system /,(z) = 2” 1s complete, we have 


Theorem 4.7. The normal system is complete with respect to every funda 


mental circle around the origin. 
We note that the coefficients a, are in this case given by the generating 


functions 


0) De U0, 2*) = ay + DF tans ieee (4.27). 
n=0 a 


4 
i 


—- 
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Table 4 
Normal Systems of Solutions of Some Differential Equations 
E ti Normal system 
berbe (apart from constant factors) 
Au =0 yt et ind 
+ind 
Au+ k?u=0 gb aers 
4 (A + 1) | LON cine 
Au 1 (i =o y2)2 0 2 (<=) é 


x=rcost, y=rsind; 
upper sign: t2,, lower sign: U2, _ 4. 


As usual we have put 


Ul, 2) =u (45= *S+) 


Z 3 214 

Since for simply connected regions the functions of the complete system of 
holomorphic functions may be taken as polynomials, the functions of a complete 
system of solutions may be assumed to be finite linear combinations of the 
normal solutions, if the latter can be defined. This means the same as 


Theorem 4.8. I} the domain © is contained in a fundamental domain which 
includes the origin, then every solution regular in © can be approximated uniformly 
in ®© by linear combinations of the normal solutions. 


4.3. The Operator Q for Singular Equations 


If a domain D contains singularities of the coefficients of the differential 
equation e(%) = O, it cannot belong to any fundamental domain. It is therefore 
not possible to define the Riemann function for such a domain, and no solution 
of e(u) = 0 which is regular in D can be generated by the method of § 4.2. For 
several differential equations with singular coefficients the operator 2 never- 
theless exists in the following sense: 

There exists a class C’ of functions /(z) holomorphic in and a class C of 
solutions u(x, y) which is contained in III (D) such that to every f(z) € C’ there 
corresponds exactly one solution in C, If u(x, y) = U(z, 2*), the function corre- 
sponding to u(x, y) is given by 


f(z) = Ul, 2) , (4.28) 


where 2* is a point in D. If we use the notation (4.17) for this correspondence, 
then the operator 2 is again linear and continuous. Hence any system which 
is complete with respect to C’ is transformed into a system complete with 
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respect to C. If z= 0 and if the complete system of holomorphic functions 
consists of a subset of the functions 2" (w = 0, 1, 2, ...), We may again speak of 
a normal system of solutions. 

In Table 5 we list some differential equations with singular coefficients for 
which the above statements have been established. Also listed are the domain 
D, the class C, 2, the functions /,,(z) of the system complete with respect to 
C’ and the corresponding solutions w,(x, y). For proofs and applications of 
some of these results, see [4], [5], [6], [7], [8], [9]. 


4.4. Addition Theorems 


In this section we suppose that the differential equation e(w) = 0, singular 
or not, admits a transformation. More precisely, we shall assume that there 
exists a class C of solutions of e(w) = 0 regular in some domain © and a map- 
ping . 
Po = 3(P) (4.29) 


[P’ = (x’, y’), P =(x, y)] of G onto itself with the following property: If 
u(P)€C, then also «'(P) = u[t(P)] € C. 

It is clear that the totality of all such transformations forms a group. We 
shall denote this group by J”. 

Let now the functions w,,(P) (7 = 0,1,...) form a complete system of solu- 
tions with respect to the class C. By this we mean that every function in C 
can be expanded in a unique way into a series of the form (4.22). Under the 
above assumptions also the functions 

mt “al (PY=4,(e(P)) (n= 07152.) (4. 30) 


n 


belong to C and can therefore be expanded into series of the same form: 
gl ECP) =>. Gin a Mog 2) ol = 0, 1...) (4.31) 


We may call these expansions ‘addition theorems’, because it turns out that 
almost all of the so-called addition theorems of the classical functions of mathe- 
matical physics are special instances of (4.31). Also .some of the equations 
listed in the Tables 4 and 5 admit transformations. These transformations 
are listed in Table 6; they all give rise to addition theorems for the corre- 
sponding complete systems. For a more detailed account of some of these we 
refer to our previous publications [4], [7]. 

On the theoretical side we mention some open problems*). By the above 
there is an addition theorem for every transformation t of J’. The coefficients 


4) The author is indebted to Professor C. LOEWNER for the substance of the following remarks. 
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Table 5 
Complete Systems of Solutions for Some Singular Equations 
e(u) = 0 D C 
Au + int oO | 
y 
|z|<a 
2 F a es 
Au + — uy+ hk? u = u(x, y) = u(x, —y) 
y 
1 : 
PPV cad 5 NS ene ez hes 
Sate Oe (L + r?)? 
|z| <a | u(x, vy) = u(x, —y) 
2v+1 Dot ge, = u(—#, y) 
jug 2htd gy 2th tg || le eVRaEl <e 
oof) cas | ux, 9) = u(e,—9) 
Z+e¢ 
i 
ik 2 1 
Aish ciple oe lik u(x, y) = u(x, —y) 
é 2 zi = u(—*, y) 
+ [40k + 2? (x? + y*)] u = 0 
x=rcos?, y=rsind, O= yi — 272 cos2B + 74, 
Table 6 
Equations Admitting Transformations 
e(u) = 0 5 C ifs 
Au + k?u=0 regular solutions translations and 


2 
Au + Stet Bu = 0 


2yvy4+1 


asin Table 5 


Au 


a Qe A 
49 


y 


rotations of the plane 


mn mM (4.31) depend on t. We introduce the infinite matrices 


A(t) = ele) Geren ae 


It is easy to see that formally the consecutive application of two elements of I” 
is mirrored in the multiplication of the corresponding matrices: 


A(at) = A(o) A(t). 


(4.32) 
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Table 5 
Complete Systems of Solutions for Some Singular Equations 


25 fin(2) u,(%, Y) (apart from constant factors) 
0 re yr C" 
r” C’ (cos 8) 
0 gn eal) A 
vr" J, lk 7) Cr (cos 8) 
0 on yn! pean oa 
a P; y "Gs “ C) (cos ®) 
0 gan ye (v, 5 
gs Ms )(cos2 8) 
(», i a1 45 ae x 2 
0. | Pe“(/1—2) PMT —23) POM + 172) PIM — 72) 
i pen (| &— 6.\0 = i+y -2 —p/2 Z 
Co) ee rer aes Ms AMA ie ake ONC 
0 ze aoa ae My yev4t/2 nlh r) 
x PB?) (G08? 2) 


e rt Tie BD) 
= a 


V+ e)?—42 2 Wes |g oe ea 


This means that from a formal point of view the matrices A(t) furnish a repre- 
sentation of the group J”. It would be of interest to discuss the conditions 
under which the product (4.32) actually exists. There arises also the question 
if the representation of J” by the matrices A(t) is faithful. 


5. The Cauchy Problem 


The Cauchy problem is usually considered in connection with hyperbolic 
differential equations. In this section we shall study some aspects of the 
Cauchy problem for the elliptic equation e(u) = 0. 


5.1. Formulation of the Problem 


Let C be a simple analytic arc in the (x, y)-plane, and let f(z) and g(z) be 
two holomorphic functions defined in a domain which contains C. Let a 
positive normal y be defined on C. We consider the problem of finding a solution 
u(x, y) of e(w) = 0 which is regular in a neighborhood of C and which on C 


satisfies the conditions 


u=f(z), =—- =8(2) (Z=x+tyEC). (5.1) 
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This problem is a typical case of the so-called Cauchy problem in the 
general theory of partial differential equations. Since the considered differ- 
ential equation is elliptic, it may be called an elliptic Cauchy problem. According 
to a wide-spread conviction, Cauchy problems for elliptic equations are not 
correctly posed. Indeed, even though under the assumption that the coefficients 
of e(w) are analytic along C there exists a unique solution of the Cauchy 
problem, it may happen that in a sense this solution does not depend conti- 
nuously on the data. This is demonstrated by the following example (essen- 
tially due to J. HADAMARD): Let C be the real axis and let e(u) = Au. We con- 
sider the sequence of data (n=0,1,...) f,(z) = (cosm z)/n, g,(z) =0. The 
corresponding solutions are easily found to be 


cosn * Cojn y (5.2) 


nN 


Un(%, y) aa 


Although /,,(z), g,(z) > 0 on-C, the sequence of solutions does not tend to zero 
in any domain of the (x, y)-plane. 

In the following sections we shall try to clarify this situation and, in 
particular, to save the continuous dependence of the solution on the data by 
subjecting the latter to stronger conditions. 


5.2. Results 


We shall need the concept of conformal symmetricity with respect to a 
given analytic arc. 


Definition 5.1. A domain © ts called conformally symmetric with respect to 
an analytic arc C, tf there exists a conformal map which transforms C into an 
interval of the real axis and © into a domain which is symmetric with respect to 
the real axis. 

The four following theorems will all be proved in section 5.3. 


Theorem 5.1. Let the fundamental domain © of e(u) be conformally sym- 
metric with respect to the arc C, and let the data f(z) and g(z) be holomorphic 
throughout ©. Then the solution of the Cauchy problem (5.1) exists and is 
regular in ©. 

This theorem admits of the following converse: 


Theorem 5.2. Let the fundamental domain © of e(u) be conformally sym- 
metric with respect to C, and let u(x, y) be a solution of e(u) = 0 which is regular 
in ®. Then the functions 


fe)=u, g)= Se Gaxtsyed) (5.3) 


are holomorphic on C and can be continued analytically into the whole of ©. 
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We may look at the solution of the Cauchy problem as an operation which 
associates with every pair of functions holomorphic in © a solution of e(w) = 0. 
We write accordingly 

u(x, y)= Off), 


where = (f(z), g(z)]. It should be noted that in this case the pair p does not 
have to satisfy a condition of admissibility. Moreover, the functions f(z) and 
g(z) are assumed to be holomorphic in the same domain ©. Apart from that the 
properties of 0 are analogous to those of 2. Theorem 5.2 shows that @ possesses 
an inverse. (An explicit formula for the inverse will be given in the proof of 
Theorem 5.2.) We also have the following analogues of the Theorems 4.4 and 4.5: 


Theorem 5.3. The operator @ is linear. 
Theorem 5.4. If /,(2) > f(z) and g,,(z) > g(z) uniformly in ©, and if p,, 
mo En) , &), then 


untformly in ®. 

The last theorem contains a sufficient condition for the continuous depend- 
ence of the solutions of an elliptic Cauchy problem on the data. In the example 
of Hadamard quoted above this condition is clearly violated, since in no domain 
of the complex plane we have uniformly 


ees 


5.3. Proof of the Theorems stated 
We split up the general Cauchy problem into the two partial problems 


(1) u=f@)> “=o 


FEN On 
Bee at tt pe), | 


We shall call these problems the first and the second Cauchy problem, and 
their solutions will be denoted by 0,[f] and @,{g], respectively. Clearly, 


O[p] = Of] + O,[¢] - (5.5) 


We first take up the case where the initial curve C consists of a segment of 
the real axis. We shall call this the reduced case of the Cauchy problem. Here 
the solution of the two partial problems can be written down explicitly in 
terms of the Riemann function. Let © be a fundamental domain which is 
symmetric with respect to the real axis, and let the functions f(z) and g(z) be 


holomorphic in G. Let, as before, 
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be the complex form of e(w) = 0. We then can construct the two functions 


Ulere* = fe (f(z) Rayz, 2, 2*) + f(2*) Ret, 2%, 22") 


Is [io {[Bu. CAS Ride, 2a (5.6) 


ae - (R(t, t, z, 2*) - R(t, t, z, zy} dt, 


Viz, a) = 45 Bi g(t) Rit, t, z, 2%) dé. (5.7) 


It is evident that the functions U and V are analytic in the cylinder domain 
(G, G). We assert that the functions u(x, y) = U(z, 2) and v(x, y) = V(z, 2) 
(zg = x +7) are solutions of the first and the second Cauchy problem, respec- 
tively. By construction the functions w and v are in class III (6). By formal 
computation, using the relations of Theorem 2.4, it is easy to show that 
E(U) = 0, E(V) = 0 and that the relations 


U (2,2) = f(z), = Uy area aaa (5. 8) 
and 
Viz, z)=0, Viz) — Volz, 2) = —2 22) (5.9) 
hold. This implies that e(w) = 0, e(v) = 0 and 
u(x, 0) = f(x), -Uol%, 0) = 0; (5.10) 
a(x, 0) = 0% wala; 0) = ea): (5.11) 


This proves the assertion and, since the uniqueness of the solution of the 
elliptic Cauchy problem is classical, Theorem 5.1 for the reduced case. The 
statements of Theorem 5.3 and Theorem.5.4 can be read from the integral 
representations (5.6) and (5.7). 

To prove Theorem 5.2 in the reduced case, we note that by VEKUA’S 
theorem there exists a function U(z,z*), analytic in (©, G), such that 
u(x, y) = U(z, 2) for z € G. We define 


f(z) = Uz, 2), g(z) =+[U,(z, 2) — U,(z, 2)) . 
The functions /(z) and g(z) are holomorphic in ©. However, for (x, 0) € G 
f(x) = u(x, 0), g(x) = a(x, 0). 


Thus, /(z) and g(z) represent the analytic continuations of u(x, 0) and u(x, 0) 
into ©. 


, 
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We now turn to the case of a general initial curve C, By assumption there 
exists a conformal transformation z = y(¢) which maps an interval J of the 
real axis of the ¢-plane onto C and a domain 6’ which is symmetric with 
respect to the real axis onto ©. The-initial conditions (5.1) can in terms of the 


function 
* > ok 
Ulz, 2*) = u(=+* fe 
soit cia ta ta es 


be written in the form 
U(9(s), 9(s)) = (p(s) (5.12) 
pe (Ui P): PO) #9) — Val $l), 9) 9) }= (9), 6.13) 
where s € J. We now substitute into E(U) = 0 new variables ¢ and ¢* defined by 
z= (0), 2* = 9(C*), (5.14) 


where (¢*) = p(¢*). The equation then transforms into an equation E’(V) =0 


for the function V(¢, ¢*) = U(g(f), y(é*)) which is of the same form and of 
which by Theorem 2.7 6’ is a fundamental domain. We write 


F(t) =f(p()), G(C) = g(@(¢)). (5.15) 


The functions F and G are holomorphic in 6’. For V(¢, €*) we get from (5.12) 
and (5.13) the boundary conditions 


V(s, s) = F(s), (5.16) 
(s,s) — Vals, s) = H(s) , (O10) 

where 
H(s) = —1 |g'(s)| G(s) (se Z). (5.18) 


We assert that the functions f(s) and H(s) can be continued analytically into 
’. This is evident for F(s). For H(s) the continuation is given by 


H(t) = —iV 9'(0) o'(0) G(0). (5.19) 


Hence V(¢, €*) must satisfy a reduced Cauchy problem. Since the data is 
continuable into 6’, the solution exists in (G’, G’). By transforming back we 
find that U(z, z*) is a solution of the original problem which is analytic in 


(G, G). This proves Theorem 5.1. . 

The proofs of the Theorems 5.3 and 5.4 follow from the corresponding 
statements for the reduced case. The transformation (5.14) preserves con- 
tinuity. 

To prove Theorem 5.2, let the solution u(x, y) be regular in the domain 6, 
which as before is assumed conformally symmetric with respect to C. By 
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VEkua’s theorem the continuation U(z, z*) into (6, G) exists. Let the function 
p(c) be the same as in the proof of Theorem 5.1, and let g~* be the inverse 
function. We introduce the function 


2 =e Diz) = og 4z)). (5.20) 


The function @(z) maps © conformally onto 6. (The point z is the ordinary 
conjugate of the conformally conjugate point of z with respect to C.) The 


inverse map is given by rs i! eae 
pk! C3 Be (Gah (5.21) 
We note the relations 
Tig\ P (p_*(2)) Gi (ay az P (P_*(2)) 5.22 
P= ga)? © O = wee) * Can 


We assert that the desired analytic continuations of the values of u and du/dy 
on C are given by 


f(z) = Ulz, 2), (5. 23) 
g(e) = i{ Uy(z, 2) VG") — U,(2, 2) VO") }. (5.24) 


Indeed, these functions are holomorphic for z€ G, since ®'(z) + 0, @'(2) + 0. 
But for z = g(s) and s € I we have z = @(s) and, in virtue of (5.22), 


'(o(3))= 2, @ (G9) = 2. 5.25 
(y(s)) , (9s) 6) (5.25) 
Hence, if the square roots are appropriately chosen, the functions /(g(s)) and 
g(y(s)) reduce to the expressions on the left of (5.12) and (5.13). This completes 
the proof of Theorem 5.2. 


5.4. Singular Cauchy Problems 


The results of section 5.2 are subjected to the condition that the considered 
domain is fundamental. As in section 4.3 for the operator Q, it has been possible 
in the case of some special differential equations to extend the above theory 
of the elliptic Cauchy problem to nonfundamental domains. We shall present 
these results here without giving proofs. Proofs of some of the results are given 
in [4] and [9]; for some related results, see [3], [11] and [13]. 

The mentioned equations are the following: 


Au+ — uy =0, | (5.26) 
Asi — tig + Hw = 0, (S727) 


2 2 
Au + - Uy + = Uy =0. (5.28) 
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Although the first equation is a special case of the latter two, it is considered 
separately, since for it the explicit results aré particularly simple. The three 
equations together will be referred to below as the special equations. 

The Cauchy problems to be dealt with in this section are singular in that the 
initial curve is a segment of the x-axis, which for y + 0 is a singular line of the 
special equations. The problems are all first Cauchy problems in the sense of 
section 5.3. It can be shown that for the initial curve y = 0 and for the above 
special equations the second Cauchy problem possesses in general no solution. 


Definition 5.2. A domain © is called suitable, if it is simply connected and 
symmetric with respect to the x-axis and if, in the case of equation (5.28), it lies 
either entirely in the halfplane x > 0 or entirely in x < 0. 

We now can state the following partial analogues of the Theorems 5.1 and 5.2: 


Theorem 5.5. Let the domain © be switable, and let f(z) be holomorphic in 
©. Then the three special equations possess for vy + —1, —2,... uniquely 
determined solutions u(x, y) which are in class III (G) and which for (x, 0) € G 
satisfy the conditions 

u(x, 0) = f(x), 2,(%, 0) =O. (5.29) 


Theorem 5.6. Let Rev > 0, let G be a suitable domain and let u(x, y) be a 
solution of a special equation which ts in class 1(©). Then the function 


f(x) = u(x, 0) 


can be continued analytically into ©. 

Combining the two theorems we also get the analogue of VEKuA’s theorem 
for suitable domains and for the special equations with Re v > 0. We do not 
know if this theorem is also valid for Rev < 0. 

Although the Riemann function cannot be used, since it is undefined for 
nonfundamental domains, it is still possible to represent the solutions of the 
singular Cauchy problem for the special equations explicitly in terms of f(z). 
For Re y > 0 these representations are written most simply in the form 


u(x, y) = [ts 9101 fx +iyt) dt. (5.30) 


The path of integration must be curved, if necessary, so that the argument of 
‘f lies always in G. The kernels K are tabulated in Table 7. 

For y = 0 the Cauchy problems are regular, so that the representation of the 
solutions in terms of the Riemann function can take place. For Rey S 0,» + 0 
the integral (5.30) is to be replaced by a contour integral of a type considered 
in [11]. 
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Table 7 
Kernels for the Representation of Solutions of Singulay Cauchy Problems 


Equation K(%#, 3 0) 
I (y+ 1/2 re 
ae 1/2 k =o wv PRs es 
ey rae Gra Dee or 
I'(v+ 1/2) 


(5.28) 


ae ; y—1 R q?(1 — 2?) 
Prag Uti aA Ha wt ge ap) 
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Zusammenfassung 


In einer Reihe von Abhandlungen hat I. N. VeEKua eine Theorie der ellip- 


tischen partiellen Differentialgleichungen mit zwei unabhangigen Variabein und 


mit 


analytischen Koeffizienten entwickelt, in welcher gewisse klassische Aspekte 
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der Theorie der holomorphen Funktionen einer komplexen Verdnderlichen ver- 
allgemeinert werden. Die Vekuasche Theorie ist. fiir den angewandten Mathe- 
matiker insofern von Interesse, als sie ihm ein in manchen Fallen wirksames 
konstruktives Verfahren zur Herstellung der Lésungen in die Hand gibt. 

Die vorliegende Arbeit ist als eime Einfiihrung in die Vekuasche Theorie 
gedacht. Es werden im ihr vor allem die Vekuaschen Resultate beztiglich der 
analytischen Fortsetzung der Lésungen und ihrer Darstellung durch holomorphe 
Funktionen zur Sprache gebracht. In einer zweiten Arbeit soll die damit eng zu- 
sammenhangende Vekuasche Theorie der Randwertprobleme dargestellt werden. 
Der Autor hat die Gelegenheit wahrgenommen, die Vekuaschen Resultate um 
einige Beispiele zu bereichern und auch einige eigene Resultate einzuflechten, 
die dem gleichen Gedankenkreis angehGren. Es sind neu oder fussen auf eigenen 
Arbeiten des Autors die Abschnitte iiber singulare Differentialgleichungen (§ 4. 3), 
Additionstheoreme (§ 4.4) und Cauchysche Anfangswertprobleme bei regularen 
und singularen Gleichungen (§ 5). 


(Received: August 22, 1956.) 


Bolzenkathode als Objekt im Elektronen-Emissionsmikroskop 
Von Ents B. Bas, Ziirich+) 


In friiheren Arbeiten ([1], [2], [3], [4]?)} wurde die vom Verfasser entwickelte 
Bolzenkathode im Zusammenhang ihrer Anwendung in Elektronenkanonen 
naher behandelt. In folgender Arbeit soll gezeigt werden, dass diese Kathoden- 
konstruktion auch im Elektronen-Emissionsmikroskop als Objekt gute Dienste 
leisten kann, insbesondere durch eine neuartige Praéparationsmethode. 

Anhand von Figur 1 wollen wir zuerst das Prinzip der Bolzenkathode mit 
einigen Worten kurz erlautern. Das als Bolzen bezeichnete Wolframstabchen 7 
ist an einem Ende in einer Klemmbacke 2 fest eingeklemmt, wahrend die 
Stirnflache am anderen Ende als Emissionsflache dient und in unserem Falle 
den zu untersuchenden Stoff 3 enthalt. Um den Bolzen herum ist konzentrisch 
eine Wolframwendel 4 angeordnet, welche von den Stromzufiihrungsstiften 5 
und 6 getragen wird. Diese Zufiihrungsstiften tragen zugleich einen Strahlungs- 
schutz, bestehend aus zwei konzentrisch zum Bolzen angeordneten Molybdan- 
zylindern 7. Kathodenblende 8 vervollstandigt die Kathode, wahrend 9 die 
erste und die zweite Blende des Immersionsobjektives darstellt. 

Die Heizung des Bolzens erfolgt so, dass die Wendel durch Stromdurchgang 
geheizt wird und als Primarkathode dient. Sie erhalt gegentiber dem Bolzen 
eine negative Spannung von einigen hundert Volt, womit die von der Wendel 
emittierten Elektronen auf den Bolzen aufprallen und ihn auf die gewiinschte 


1) Institut fiir technische Physik der ETH. = 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 212. 
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Temperatur aufheizen. In der Heizschaltung ist dafiir Sorge zu tragen, dass 
infolge des Riickheizeffektes von Bolzen auf Wendel keine Heizinstabilitat 
auftritt. Bei niedrigeren Objekttemperaturen kann die Bombardierung vollig 
ausbleiben. Die Heizung des Bolzens erfolgt dann durch die Zustrahlung der 
Heizwendel. 


Figur 1 


Schnitt durch die Bolzenkathode als Objekt im Emissionsmikroskop: 7 der Bolzen aus Wolfram als 
Objekthalter; 2 Einspannbacken fiir den Bolzen; 3 auf den Bolzen aufgeschmolzenes Objekt; 
4 Heizwendel aus Wolfram; 5 und 6 Haltestiften fiir die Halterung der Heizwendel und der Ab- 
schirmung; 7 Abschirmzylinder aus Molybdan; 8 Abdeckkappe aus Molybdan; 9 Lochblenden des 
Immersionsobjektivs. 
Die Heizung des Objektes erfolgt so, dass die Heizwendel 4 mit Gleichstrom auf etwa 2300°C 
geheizt wird; sie emittiert dann Elektronen, welche auf rund 300 V beschleunigt werden und den 
Bolzen 1 bombardieren, wodurch seine Temperatur und somit die Temperatur des Objektes 3 auf 
den gewtinschten Wert erhoht wird. 


Fir die Anwendung der Bolzenkathode als Objekt im Emissionsmikroskop 
sind vor allem drei Dinge von Bedeutung: Axiale Ausdehnung des Bolzens, 
Schwingungen des Bolzens und das magnetische Stérfeld der Heizwendel. 

Die axiale Ausdehnung des Bolzens ist nicht zu vermeiden. Dies hat zur 
Folge, dass man bei Veranderung der Objekttemperatur die Scharfstellung 
durch Verschiebung des Immersionsobjektives nachkorrigieren muss. Dieser 
Umstand ist besonders bei raschen Temperaturanderungen, welche durch die 
Dynamik der zu untersuchenden Vorgaénge vorgegeben werden, stérend. Diese 
thermische Bewegung des Objektes ist aber bei allen Objektanordnungen in 
den Emissionsmikroskopen mehr oder weniger vorhanden. Die Emissionsflache 
der Bolzenkathode bei der in Figur 1 dargestellten Konstruktion verschiebt 
sich beim Aufheizen auf 2000° K Objekttemperatur um rund 0,1 mm, und bei 
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1000° K betragt die Verschiebung etwa 0,04 mm. Die Verschiebung des Objektes 
bei konstanter Objekttemperatur infolge Nichterreichen des Beharrungszu- 
standes ist meistens so gering, dass Belichtungen bis zu 30s Belichtungszeit 
dadurch nicht gestért werden. 

Den mechanischen Schwingungen des Bolzens ist grosse Bedeutung beizu- 
messen, da sie die erzielbare Auflésung stark beeinflussen kénnen. Um hier 
einen Uberblick zu bekommen, betrachten wir den idealisierten Fall eines ein- 
seitig eingespannten elastischen Stabes von der Linge / mit einem Massen- 
punkt m am Ende (Figur 2). Die Transversalkoordinate dieses Massenpunktes 


Figur 2 


Mechanische Schwingungen des Bolzens. 


bezeichnen wir mit y* und die laufenden Koordinaten des Stabes mit x, y. 
Dann lasst sich die Schwingungsgleichung dieses Systems wie folgt formulieren : 
d*y* 


ica 


ad? h(a? 
(— 4) FEIZA +81 = (G4) =9- (1) 
Das letzte Glied, das der zeitlichen Anderung des Biegemomentes, also des 
zweiten Gliedes, proportional gesetzt worden ist, beriicksichtigt die Dampfung 
durch innere Reibung. & ist der sogenannte normale Viskositatskoeffizient [5]. 
I bedeutet das equatoriale Trigheitsmoment des Stabes und £ den Elastizitats- 
modul. Die zweimalige Integration obiger Gleichung nach x mit Einsetzen der 


Grenzbedingungen x = 0; dy/dx = 0, y=0O und x =1; y= y* liefert uns: 


d*y 


<2 dy* 3 Pe 
m—a + a lt gt Ly = 0), (2) 


Dies ist eine Gleichung fiir gedampfte Schwingungen, deren Lésung lautet: 


yt = yt e Ft? cos (am, t+ e) (3) 
mit Ec 
nw es 

Wo = Serie : (4) 


als Eigenfrequenz des Systems. | 
Fiir die Halbwertzeit der gedampften Schwingungen erhalten wir: 


lens ge (5) 


206 Ents B. Bas ZAMP 


Fiir die Anwendung dieser Betrachtungen auf die Bolzenkathode wollen 
wir folgende Vereinbarung treffen: als die Stablange / soll derjenige Abschnitt 
gelten, welcher sich nicht auf hoher Temperatur befindet und dadurch seine 
Elastizitat bewahrt. Als Masse m gilt die Masse des tibrigbleibenden Endstiickes 
des Bolzens. Nehmen wir als Beispiel fiir unsere in Figur 1 dargestellte Kon- 
struktion mit: 


t= 7-+10-'cm; dz,=6-10cem; 1,—3= 60°10 tem: 4m = 32/7 eee 
und fiir Wolfram 
E = 3,5 «101 Dyn/em?;, 6 10° Dyn eam; 


so erhalten wir: 
Oy = 22,9 kHz; iw 10-* s. 


Nehmen wir weiterhin an, dass der Bolzen durch einen Stoss auf eine Amplitude 
von: 
ye = 10-7 cm 


angeregt wiirde, so wird diese Amplitude in 10-3 s auf 10-2. 10-3 = 10-®° cm 
oder 0,1 w absinken. Wir sehen hieraus, dass die Erschiitterungsempfindlichkeit 
bis herab zu der Aufldsungsgrenze des Emissionsmikroskopes ausreichen diirfte. 

Die dritte uns hier interessierende Frage betrifft das magnetische Stérfeld 
im Immersionsobjektiv, hervorgerufen durch den Wendelheizstrom. Das magne- 
tische Feld der Wendel kénnen wir in zwei Anteile zerlegen, den rotationssym- 
metrischen Anteil, hervorgerufen durch die Windungen der Wendel, und den 
Querfeldanteil, hervorgerufen durch die Wendelzuleitungen. 

Die rotationssymmetrische Komponente ist sofort aus der Feldstarkenglei- 
chung einer Zylinderspule anzugeben: 


eee. 


| Hatneliuigs achull = eee 
” r+) J G=y | 


I,, Wendelheizstrom; 

w Anzahl der Wendelwindungen; 

l, Wendellange; 

d, mittlerer Wendeldurchmesser ; 

z laufende Koordinate auf der Achse der Wendel, von der Mitte der Wendel 
gezahlt. 


ape 


Hier bedeutet: 
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In unserer Anordnung hat die Heizwendel folgende Daten: 


lo = 3,5 mr 
@g=1,35 mm; 


(aaa 

Pes 2A 

damit erhalten wir fiir die axiale Feldstarke die Beziehung: 
et ee ae ee [Oe] 
PeCRES Ah\2 iG Sheer LISS 
1 CMS ie deat picked sek 
|/ + (spas) eee ae 
[zin mm]. 


Diese Beziehung ist in Figur 3 graphisch dargestellt. 


"4 
2 
2 
o 
A=} 
10 
> 
eirb 
a 
a Iy= 2A 
6 dw=35mm 
w =14 
lh = 4mm 
4 
2 
0 


fie Dees Olle Game O ozo 
Figur 3 
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Axialsymmetrische und transversale Anteile der magnetischen Stérfeldstarke der Heizwendel. 


Die Querkomponente des magnetischen Wendelstérfeldes berechnen wir 
aus dem in Figur 4 dargestellten vereinfachten Stromkreis. Nach dem Biot- 
Sovartschen Gesetz lasst sich fiir jede einzelne Stromkreisstrecke fiir den in 


einem Punkt P erzeugten Feldanteil angeben: 


ay 
I sin « 
SP alte cin mies 
oy 


(8) 
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Figur 4 


Zur Berechnung der axialen und radialen magnetischen Stérfeldstarke der Heizwendel. 


In unserem Falle liegen alle Leiter in der Zeichnungsebene; damit ist die Feld- 
starkenrichtung stets zur Zeichenebene senkrecht. Da wir uns nur fiir die Feld- 
starke auf der Achse interessieren, liefern nur die Strecken 7, 2, 4 und 3 einen 
Beitrag. Die Rechnung fiihrt schlussendlich auf folgende Beziehung: 


B= a WVi+ Cae) ++ Gael -2] Wem. 


Setzen wir wieder unsere obigen Wendeldaten und fiir die Lange des Zuleitungs- 
drahtes J, = 4mm ein, so erhalten wir: 


Ns V2 + ( 5 pe 5) pee le x (a) = 2| [Oe] (10) 


[gin mm]. 


In Figur 3 finden wir auch H, in Abhangigkeit von z graphisch aufgetragen. 
Aus konstruktiven Griinden betragt der minimale Abstand der Emissionsflache 
von der Wendelmitte ungefahr 2,8 mm. Fiir diesen Fall erhalten wir fiir die 
beiden Feldstarkenkomponenten folgende Werte 


H,=7,3(Ocly H, = 1:63(Oe). 


Der Einfluss dieses Stérfeldes auf die Abbildungsgiite des Immersionsob- 
jektives ist schwer zu beurteilen. Fallt die Achse der Heizwendel mit der Achse 
des Immersionsobjektives zusammen, so ditirfte die Wirkung der rotationssym- 
metrischen Feldkomponente H, ganzlich zu vernachlassigen sein. Hingegen ver- 
ursacht die transversale Komponente H, stets eine ablenkende Wirkung auf 
die emittierten Elektronen, so dass die optische Achse gekriimmt wird, wodurch 
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zusatzliche Abbildungsfehler entstehen kénnen. Selbstverstiandlich wird die 
Wendel mit gutgesiebtem Gleichstrom geheizt, so dass keine Wechselkompo- 
nente des Feldes vorhanden ist. Die obige Angabe fiir den Abstand zwischen 
Wendelmitte und Emissionsflache von 2,8 mm ist der konstruktiven Ausfiih- 
rung der Bolzenkathode als Elektronenquelle fiir die Elektronenkanonen ent- 
nommen. Bei dieser Anwendung ist man bestrebt, die Bombardierungszone 
méglichst nahe bei der Emissionsflache anzuordnen, da sich die heisseste Stelle 
des Bolzens in der Mitte der Bombardierungszone befindet und durch die 
starke Verdampfung des Wolframs in dieser Zone die Lebensdauer der Kathode 
bedingt ist [4]. Nun sind wir beim Emissionsmikroskop aus zwei Griinden in 
dieser Beziehung nicht so stark eingeschrankt. Erstens ist die Lebensdauer des 
Bolzens von untergeordneter Bedeutung, das heisst, wir kénnen fiir die Bolzen- 
temperatur in der Mitte der Bombardierungszone viel héhere Werte zulassen, 
und zweitens sind Temperaturen der Emissionsflache bei den elektronenmikro- 
skopischen Untersuchungen bedeutend tiefer als bei der Anwendung der 
Bolzenkathode als Elektronenquelle in Elektronenkanonen. Aus diesem Grunde 
kann die Emissionsflache bedeutend weiter von der Heizwendel weggeriickt 
werden, wodurch nach Figur.3 die Stérfeldstarke der Heizwendel noch weit- 
gehend vermindert werden kann. Ausserdem erméglicht ein grosser Abstand 
zwischen Heizwendel und Emissionsflache die Anordnung einer magnetischen 
Abschirmung. Es erhebt sich nun die Frage, wie sich die Bombardierungs- 
leistung des Bolzens fiir eine bestimmte Objekttemperatur beim Wegriicken 
der Emissionsflache von der Heizwendel verandert. Diese Frage wurde experi- 
mentell gepriift, woraus die Kurvenschar in Figur 5 resultiert. Die Messungen 


Parameter: Temperatur der Emissionsflache 


5 


10 


Bombardierungsleistung \w) 


25 3,0 35 40 45 5,0 55 
Abstand Heizwendelmitte— Emissionsflache [mm] 


Figur 5 


Abhangigkeit der Heizleistung von der Entfernung der Emissionsflaiche von der Heizwendelmitte: 
Bolzendurchmesser 0,6 mm (Wolfram), Objektdurchmesser !,5 mm (Wolfram). 


ZAMP VIII/14 
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wurden mit einem W-Objekt durchgefiihrt, und die Anordnung war so getroffen, 
dass der Abstand zwischen Wendelmitte und Einspannstelle des Bolzens stets 
11,2 mm betrug. Bolzendurchmesser war dz = 0,6 mm und Objektdurchmesser 
dx = 1,5 mm. Durch die Verlangerung des hochgeheizten Bolzenendes wird nun 
auch die axiale Ausdehnung des Bolzens ausgepragter. Die gemessenen Aus- 
dehnungen sind durch die Kurvenschar in Figur 6 graphisch dargestellt. Uber 
den Einfluss des Abstandes Wendelmitte-Emissionsflache auf das Auflésungs- 
vermégen wurden noch keine systematischen Versuche durchgefuhrt. 


200 
Parameter: Temperatur der Emissionsflache 


tu] 
Fite 


Ausdehnung des Bolzens 


25 3,0 35 40 45 $,0 5,5 
Abstand Heizwendelmitte— Emissionsflache [mm 


Figur 6 


Axiale Verschiebung der Emissionsflache bei verschiedenen Temperaturen in Abhangigkeit des 
Abstandes Emissionsflache—Heizwendelmitte: Bolzendurchmesser 0,6 mm (Wolfram), Objektdurch- 
, messer 1,5 mm (Wolfram). 


Wir wenden uns jetzt der Objektpraparation zu. Diese besteht aus drei 
Stufen: 
a) Praparation des Wolframbolzens; 
b) Aufschmelzen des zu untersuchenden Stoffes auf den Wolframbolzen; 
c) Praparation der Emissionsflache. . 


Der Wolframbolzen wird aus dem handelsiiblichen Wolframdraht herge- 
stellt. Der meist angewendete Durchmesser betragt 0,6 mm. Zuerst werden 
Drahtstticke von rund 300 mm Lange in einer Streckvorrichtung in Schutzgas 
bei etwa 1000°C gestreckt. Aus so gerichtetem Draht werden Stiicke von pas- 
sender Lange herausgeschnitten und in einem Rekristallisationsofen bei 3000° K 
rekristallisiert (Figur 7a). Nun werden die rekristallisierten Bolzen in einen 
Halter eingesteckt und in eine Kathodenstrahlréhre eingeschleust. Durch den 
auf die Stirnflache des Wolframbolzens konzentrierten Elektronenstrahl wird 
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das Wolfram so weit geschmolzen, bis eine“Kugel von gewiinschtem Durch- 
messer entsteht (Figur 7b), Es stand uns im Elektronenstrahl eine Leistung von 
14kV-10 mA = 140 W zur Verfiigung, womit Kugeldurchmesser von etwa 
2mm ohne Schwierigkeit erzielt werden kénnen. Schliesslich wird die Wolfram- 
kugel bis auf die Halfte abgeschliffen, und der Wolframbolzen steht fiir weitere 
Praparation des Objektes bereit (Figur 7c). 


Figur 7 


Verschiedene Stadien der Objektpraparation: a Ausgangsstabchen fiir Wolframbolzen; 6 Wolfram- 

stabchen mit geschmolzener Kugel an der Stirnfliche; c fertig praparierter Wolframbolzen; 

d Bolzen_ mit einem aufgeklebten Praparat; e Bolzen mit einem aufgepressten Pulverpraparat; 

f Bolzen mit angeschmolzenem Praparat; g einsatzbereites Objekt mit geschliffener und feinpo- 
lierter Emissionsflache. 


Der nachste Schritt ist das Aufschmelzen des zu untersuchenden Stoffes. 
Dieser Stoff kann in massiver Form oder in Pulverform vorliegen. Liegt er in 
massiver Form vor, so wird ein Stiick von geeigneter Grosse auf die Stirnflache 
eines wie oben vorbereiteten Wolframbolzens mit Hilfe von Kollodium ange- 
klebt (Figur 7d). Falls der zu untersuchende Stoff in Pulverform vorliegt, wird 
eine bestimmte Menge von diesem Pulver bzw. Pulvergemisch in einer kleinen 
Pressform auf die Planflache der Halbkugelkalotte des Bolzens aufgepresst 
(Figur 7e). Es kann ein wenig Kollodium als Bindemittel zugegeben werden. 

Der Wolframbolzen mit aufgeklebtem bzw. aufgepresstem Praparat wird 
wieder in den Elektronenschmelzofen eingeschleusst und durch Elektronen- 
beschuss auf das Wolfram aufgeschmolzen (Figur 7/). Da Wolfram einen hohen 
Schmelzpunkt besitzt und das Aufschmelzen unter Umstanden in Bruchteilen 
einer Sekunde erfolgt, besteht sehr geringe Wahrscheinlichkeit fiir eine Legie- 
rung mit Wolfram. 

Schliesslich kommt die Fertigstellung des Objektes. Dafiir wird zuerst die 
Schmelzkugel auf dem Wolframbolzen auf den gewiinschten Durchmesser 
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rundgeschliffen und dann die Emissionsflache plangeschliffen (Figur 7g). Bei 
diesem letzten Schritt ist zu beachten, dass die resultierende Dicke des zu unter- 
suchenden Stoffes auf dem Wolfram ausreichend ist. Die Schlussphase der 
Praparation bildet die feine Politur der Emissionsflache. 

Die beschriebene Objektpraparationsmethode bietet verschiedene Vorteile. 
Das Objekt kann sehr prazise hergestellt werden. Das Schmelzen erfolgt in 
kurzer Zeit, so dass keine Reaktion mit der Unterlage, das heisst mit dem 
W-Bolzen zu befiirchten ist. Man kann auf diese Weise aus kleinen Stoffmengen 
sehr reine Praparate herstellen. Dass man dabei von Pulver ausgehen kann, 
ist ein weiterer Vorteil. Man kann leicht legierte Objekte herstellen und dabei 
auch Legierungskomponenten mit héherem Dampfdruck anwenden, da der 
Schmelzvorgang sehr rasch erfolgt. Beziiglich der Schmelztemperatur ist man 
nicht beschrankt und kann mit dem Schmelzpunkt des zu untersuchenden 
Stoffes bis nahe an die Schmelztemperatur des Wolframs herangehen. Einen 
gewissen Nachteil bedeutet die Einschrankung im Objektdurchmesser. Objekte 
mit einem Durchmesser der Emissionsflache von 2 mm kénnen ohne Schwierig- 
keit hergestellt werden. 

Selbstverstandlich kénnen in vielen Fallen Bolzen aus vollem Material des 
zu untersuchenden Stoffes herausgearbeitet und in die Bolzenkathoden- 
Anordnung als Objekt eingebaut werden. Wir haben allerdings zurzeit diesbe- 
zligliche Erfahrungen nur mit hochschmelzenden Metallen Molybdaén und 
Tantal gemacht. 

Zum Schluss zeigt uns Figur 8 als Beispiel ein emissionsmikroskopisches 
Bild, wobei die beschriebene Anordnung mit einem Molybdan-Objekt zur An- 
wendung kam. Die Aufnahme stellt die Auskristallisation von Mo,C in mit C 
iibersattigtem Mo dar. Uber diese Untersuchungen wird an anderer Stelle aus- 
fiihrlich berichtet [6]. 


LITERATURVERZEICHNIS 


[1] E. B. Bas, Mitteilungen des Instituts fiir technische Physik an der Bee 
sischen Technischen Hochschule Ziirich, Nr. 1 (1950). 

[2] E. B. Bas, G. F. F. Mitteilungen 70, 17 (1954). 

[3] E. B. Bag, Optik 72, 377 (1955). 

[4] E. B. Bas, Z. angew. Phys. 7, 337 (1955). 

[5] R. Honpa und S. Konno, Phil. Mag. 42, 121 (1921). 

[6] E. B. Bas, W. Epprecut und L. Preuss, Acta Metallurgica. 


Summary 


Experiences are reported in connection with the use of bolt cathode as the 
object-holder in an electron emission microscope. The bolt cathode consists of a 
cylindrical Tungsten rod, the bolt, close pinched at one end, the top of the other 
free end being used as the emission surface. The bolt is surrounded by a Tungsten 
helix, which can be heated by passing direct current. This helix serves as a 
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1800°K 


Figur 8 
Auskristallisation von Mo,C aus mit Kohleustoff gesittigtem Molybdan. 


primary cathode emitting electrons for electron bombardment of the bolt. In 
the case of use of bolt cathode as specimen holder in an electron emission micro- 
scope the specimen is melted on the top of the Tungsten bolt. Melting occurs by 
the bombardment of specimen substances with accelerated electrons from an 
electron gun. This method has many advantages and is described in detail. Two 
properties of the bolt cathode, important for the emission microscope, i. e. the 
axial displacement of the object and the magnetic field of the heating helix, are 
discussed in detail. A micrograph shows the quality of the electron emission 


pictures obtained. 


(Eingegangen: 6. Dezember 1956.) 
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Das Verfahren det Treppeniteration und verwandte Verfahren 
zut Losung algebraischer Eigenwertprobleme 


Von FriEpRICcH L. BAUER, Miinchen?) 


In Verallgemeinerung des klassischen Bernoullischen Verfahrens haben wir 
vor kurzem Verfahren angegeben [2]?), die unmittelbar zu einer Faktorisierung 
eines vorgegebenen Polynoms P(x) fiihren. Diese Verfahren lassen in Matrix- 
schreibweise [unter Verwendung der zu P(x) gehérigen Frobenius-Matrix] un- 
schwer eine Verallgemeinerung auf den Fall eines Eigenwertproblems fiir eine 
beliebige Matrix A erkennen. Die Iterationsmatrix J’, aus / Spalten ist von 
«Trapezgestalt», I’,,, entsteht aus A J’, durch deren Orthogonalisierung an 
den ersten / Zeilen der Einheitsmatrix («Treppeniteration»). Es zeigt sich, dass 
diese Massnahme bereits zur Stabilitat fiihrt. 

Die erwahnte Orthogonalisierung (auch an einem allgemeineren System C 
von festen Zeilenvektoren, «stabilisierte Simultan-Iteration») fihrt auf die 
Dreieckszerlegung von G; = CA*B, wobei B die Ausgangsnaherung ist. Dabei 
treten Quotienten gewisser Minoren von G; auf. Solche Ausdriicke haben bereits 
PERRON [6] und spater MUntz [5] sowie AITKEN [1] studiert. Es zeigt sich aber, 
dass die stabilisierte Iteration und insbesondere ihr praktisch wichtigster Spe- 
zialfall, die «Treppeniteration», diese Minoren rekursiv und ohne Ausléschung 
liefert, im Gegensatz zu ihrer direkten Berechnung aus G;. Man kénnte die 
Treppeniteration insofern als eine numerisch vorteilhafte Abwandlung der 
Perronschen Minorenmethode bezeichnen, mit der sie theoretisch parallel layft. 
Die (in ihren Bedingungen tibersichtliche) Konvergenztheorie lasst sich weit- 
gehend auf die erwahnten alteren Untersuchungen zuriickfiihren®). 

Es liegt nahe, die Iteration an Spalten wnd an Zeilenvektoren durchzu- 
fiihren und diese gegenseitig zu orthogonalisieren (und damit zu stabilisieren). 
Die darauf aufbauende Methode der «gegenseitig orthogonalisierten simultanen 
Bi-Iteration» hat merkwiirdigerweise enge Verwandtschaft zum erstgeschilder- 
ten Verfahren. Man kénnte sie aber auch als vorschnellende Abart der Koch- 
schen Stabilisierung [4] deuten. 

Beide Verfahren stehen in engen theoretischen Beziehungen zu RUuTIS- 
HAUSERS LR-Transformation [9], zur A-Transformation von RUTISHAUSER 
und dem Verfasser [7] und zu einer kiirzlich vorgeschlagenen Eigenvektor- 


1) Mathematisches Institut der Technischen Hochschule Miinchen und Arbeitsgruppe fiir 
elektronische Rechenanlagen der Technischen Hochschule Miinchen. 

*) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 234. 

3) Vgl. Anhang, A 2, Anmerkung 22), 
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bestimmung mittels der Z R-Transformation [3], All diese Verfahren beruhen 
auf dem Grenzverhalten gewisser, einer Bernoullischen Differenzengleichung 
gehorchender Folgen («Bernoullischer Konvergenztyp») [3]. 

Die Treppeniteration und die Bi-Iteration konvergieren, abgesehen von 
singularen Fallen, unabhangig von spezieller Wahl von B und C, erfordern also 
keine sogenannte «gute Ausgangsnaherung». Sie sind stabil und selbstkorrigie- 
rend; der Rechengang ist weitgehend zyklisch (Rechenvorschrift A, Rechen- 
vorschrift B). Die Verfahren diirften damit fiir programmgesteuerte Rechen- 
automaten besonders geeignet erscheinen. 


§ 1. Treppeniteration 
1.1. Grundlegende Beziehungen 


A sei eine ” x m-reihige Matrix, von der / Eigenwerte (1 </ <n) und even- 
tuell die zugehdrigen Zeilen- oder Spalteneigenvektoren zu bestimmen seien. 
Entsprechend wird die Iteration an / n-dimensionalen Zeilen- oder Spalten- 
vektoren, die jeweils zu einer rechteckigen Matrix J", (u = 0,1, 2, ...) zusam- 
mengefasst seien, geschehen. Wir beschranken uns in der expliziten Schreib- 
weise zunachst auf den Fall der Spalteniteration, aus dem durch Transposition 
der Fall der Zeileniteration hervorgeht. 

Um zu erzwingen, dass die fn im Laufe der Iteration stets einen /-dimensio- 
nalen Raum aufspannen, verlangen wir, dass oberhalb der Hauptdiagonale 
die Elemente eines jeden J’; verschwinden*). Wir fiihren ferner die Normierung 
ein, dass die Diagonalelemente von J’; Einsen sind. 

I’; sei von dieser Trapezgestalt. [’;,, soll dann dadurch entstehen, dass 
A I’; durch geeignete Linearkombinationen der Spalten in die verlangte Gestalt 
gebracht wird, das heisst durch Elimination unter Abspaltung einer / x /-rei- 
higen oberen Dreiecksmatrix FR; ,,: 


AD; = Ti41 Riga: (1) 


Wir wollen im Augenblick offenlassen, ob solche R, stets existieren, k6nnen 
aber bereits sagen, dass / nicht grésser sein darf als der Rang m von A”, anderen- 
falls die Trapezgestalt, das heisst / linear unabhangige Spalten, sicher nicht 
standig wiederhergestellt werden kann (letzteres wiirde sich aber bereits zu 
Beginn der Iteration bemerkbar machen). Wir kénnen ferner annehmen, dass 
bereits die Startmatrix J’, der Iteration durch Zerlegung einer » x /-reihigen 
Ausgangsmatrix B gewonnen wird: 


B=TI,Ry. (2) 


4) Anders ausgedriickt, verlangen wir, dass die k-te Spalte von I’ , orthogonal zu den ersten k—1 
Zeilen der Einheitsmatrix ist. Andere Orthogonalitatsforderungen werdeu wir in § 3 betrachten. 
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Wir werden die Iteration als konvergent bezeichnen, wenn J’, fiir 4 >oo 
elementweise konvergiert und damit einen Grenzwert besitzt: 


C,>I (u>oo). (3a) 


Unter welchen Bedingungen Konvergenz herrscht, werden wir weiter unten 
diskutieren. Falls die I’, konvergieren, konvergieren aber auch die R,,, es gilt 
dann 

Wile ale (alm) (3b) 
und 


Ar=TR. (4) 


Das bedeutet, dass J’ einen Teilraum aufspannt, der durch A in sich transfor- 
miert wird. Die Eigenwerte von R sind also Eigenwerte von A. Da jedoch R 
eine (obere) Dreiecksmatrix ist, sind seine Diagonalelemente charakteristische 
Zahlen, im Falle linearer Elementarteiler sind es / Eigenwerte von A. 

Die (Spalten-)Eigenvektoren x, von A ergeben sich zu 


i=l pe (eel cases (5) 


wo 1, Eigenvektoren von R sind. Entsprechendes gilt fiir Hauptvektoren. 

Da R von Dreiecksgestalt ist, bietet die Bestimmung der n, durch Herauf- 
rechnen kein weiteres Problem. Die Matrix { 9, } (o=1,2,..°, 1) der Eigen- und 
Hauptvektoren von R ist eine obere Dreiecksmatrix. Demnach ist die Matrix 
{ x, } (o=1,2,..., 2) vonl Eigen- (oder Haupt-) Vektoren von A gegeben durch 


{z,}=I{»,} (6 


als Produkt einer unteren Trapez- und einer oberen Dreiecksmatrix. Wir 
kénnen also sagen: Die Treppeniteration liefert im Konvergenzfall / Eigen- (und 
Haupt-) Vektoren in Dreieckszerlegung. 

Fir /=1 geht die Treppeniteration trivialerweise in die Vektoriteration 
[6], [5] tiber. 


1.2. Rechenvorschrift der Treppeniteration 


Wir fassen die Rechenvorschrift in algorithmischer Schreibweise5) zu- 
sammen: 
Rechenvorschrift A. Gegeben A und die Ausgangsmatrix B= B, von 1 


Spalten: 
Fiir: = 0(1) N: triang(B) >T;,, 


ATA Seas 


5) Das gerichtete Gleichheitszeichen deutet das numerische Ergebnis an. Weitere Erklarungen 
der Symbole siehe unten. 
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t=N, falls || T;— Fy] Se, 


trrang( Bye) Pray Rua > 
diag(Ry ,,) >{é}. 

rep Ap) > Do » 

Py 410¢ > ha 


Der mit triang(B,) bezeichnete Prozess kann entweder durch GauBsche 
Elimination erfolgen, da R; zunachst nicht gebraucht wird, und schliesslich die 
Diagonalelemente von Ry, [mit diag(Ry,,) bezeichnet] die Reziproka der 
Eliminationselemente sind. Er kann aber auch durch explizite Dreieckszerle- 
gung nach BANACHIEWICZ erfolgen, wobei die R; nicht gespeichert zu werden 
brauchen. 


| C || bezeichnet eine geeignete Norm, etwa Max(c in). Das Abbrechkrite- 


Fiirea@=1(1)/: eigenv(R, 


rium kann verscharft werden, etwa dahingehend, ae die Iteration tiber einige 
Schritte hinweg «stehen» soll. Die Bestimmung des Eigenvektors durch Herauf- 
rechnen an der Matrix Ry ,, — A, E ist mit eigenv(Ry ,1, A.) bezeichnet. 


§ 2. Konvergenz 


2.1. Vorbemerkungen 


Es mag angebracht sein, die Konvergenzuntersuchungen in den Anhang zu 
verweisen. Um eine richtige Einschatzung der Konvergenzbedingungen zu er- 
méglichen, muss jedoch einer Zusammenfassung der Ergebnisse ein Uberblick 
uber das Funktionieren der Konvergenz vorangestellt werden. 

Nach (1), (2) gilt 

AP Bes Fy Pe; (7) 
wo 
Py Re Reyes Ry Ry (8) 


ist. I’, ist also die linke Dreieckszerlegte von A* B. Diese Beziebung, die die 
Grundlage der allgemeinen Konvergenzuntersuchungen von A 2 ist, lasst 
bereits erkennen, dass bei der Treppeniteration Minoren von Matrixpotenzen 
eine Rolle spielen*). Das bringt mit sich, dass sie die Tendenz hat, mit linearer 
Konvergenz eine Trennung der Eigen- (und Haupt-) Vektoren herbeizufiihren 
derart, dass diese sich geordnet nach Betragen der zugehorigen charakteristischen 
Zahlen ergeben. Das heisst, dass I"; gegen die untere Dreiecks-(Trapez-)Zerlegte 
I’ der in eben dieser Reihenfolge aneinandergereihten Eigen- und Haupt- 
vektoren streben méchte. 


6) Vergleiche die historischen Bemerkungen in der Einleitung. 
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Die oben erwahnte Ordnung gelingt selbstverstandlich nur insoweit, als die 
betreffenden charakteristischen Zahlen betragsverschieden sind. Jedoch besteht 
auch noch in dem wichtigen Fall, dass die betragsgleichen charakteristischen 
Zahlen sogar gleich sind, eine gewisse Tendenz zur Trennung von zugehorigen 
Eigenvektoren’). Diese normale Konvergenztendenz kann durch besondere Wahl 
von B gestért werden®); die Falle anomaler Konvergenz, in denen sich die 
Eigen- (und Haupt-) Vektoren nicht in der normalen Ordnung ergeben, sind 
durch das Verschwinden der Ausdriicke (*), A 2, gekennzeichnet. Praktisch 
interessiert nicht die explizite Gestalt dieser Bedingungen, sondern nur der 
Umstand, dass es sich dabei um isolierte Falle handelt®). Wenn nicht besondere 
Umstinde vorliegen?®), ist ihr Auftreten zufallig und praktisch bedeutungslos, 
selbst unter Beriicksichtigung des Umstands, dass bereits die Umgebung 
solcherart singularer B zu den erwahnten numerischen Unzutraglichkeiten 
fiihrt, also bei zufalliger Wahl von B eine endliche, aber geringe Wahrschein- 
lichkeit fiir die praktische Stérung der normalen Konvergenztendenz besteht. 

Unter der normalen Konvergenztendenz jedoch ist die Iteration numerisch 
stabil und in Verbindung mit der Tatsache, dass die Rechenvorschrift die Form 
einer Stationaritatsbedingung hat, selbstkorrigierend in der Gewinnung von /™?). 

Eine andere Stérung der numerischen Rechnung kann dadurch hervorge- 
rufen werden, dass, anschaulich gesprochen, die Dreieckszerlegung der Matrix 
Lie von Eigen- (und Haupt-) Vektoren, die sich unter einer bestimmten Kon- 
vergenztendenz, insbesondere der normalen, einstellen méchten, nicht existiert. 
In diesem isolierten Fall, der wiederum durch das Verschwinden gewisser Aus- 
driicke (**), vgl. A 2, gekennzeichnet ist und praktisch ebensowenig Bedeutung 
besitzt wie die Falle anomaler Konvergenztendenz, liegt eine ungliickliche 
Orientierung der Eigenvektoren im Hinblick auf die starre Eliminationsreihen- 
folge vor. Er kann durch Pivotelimination stets umgangen werden (vgl. 3.2), 
eine allgemeinere Behebung dieser Schwierigkeit wird sich in § 3 ergeben. 


2.2. Konvergenzsdtze 


Alles fiir die praktische Anwendung Wesentliche ergibt sich also aus dem 
nachfolgenden Satz, bei dem, wie im folgenden, die charakteristischen Zahlen 


*) Von Hauptvektoren ebenfalls, jedoch ist die Konvergenz dann nur noch logarithmisch und 
damit praktisch bedeutungslos. Es mag nur vermerkt werden, dass sich zuerst die Eigenvektoren, 
dann die Hauptvektoren, geordnet nach ihrer Stufe, einstellen. 

8) Ist etwa jede der J Spalten von B bereits Eigenvektor, so ist die Treppeniteration trivialer- 
weise sofort stationar und liefert eben diese Eigenvektoren in der urspriinglichen Reihenfolge. 

8) Sie fihren zu Unannehmlichkeiten in der numerischen Rechnung durch das Auftreten von 
Stellenverlust infolge Ausléschung; sie sind auch numerisch, das heisst unter der Auswirkung von 
Rundungsfehlern, instabil. Die normale Konvergenztendenz setzt sich also schliesslich immer durch. 
Praktisch ist dies nicht sehr bedeutungsvoll, da es unertraglich lange dauern kann. 

10) In Fallen absichtlicher Herbeifiihrung anomaler Konvergenz ist jedoch auf Stabilisierung 
zu achten °), 

11) Die schliessliche Berechnung der Eigenvektoren mag jedoch unerwiinschte Rundungsein- 
flusse mit sich bringen. Abhilfe ergibt das Verfahren von § 4. 


Vol. VIII, 1957 Verfahren zur Lésung algebraischer Eigenwertprobleme 219 
von A als betragsmissig geordnet anzunehmen sind: 


Satz 1. Die I’;, Rechenvorschrift A, konvergieren (abgesehen von singuldren 
Fallen), jalls fir k =A, 2, ...,1 


|E;.| =. lEx+1| oder §,= re 


gilt. Die Bedingungen entfallen fiir k = n. Die &, ... &, treten in dieser Reihenfolge 
mm der Diagonale von R auf. Der Iterationsprozess ist numerisch stabil und selbst- 
korrigterend. Die Konvergenz ist linear, falls alle mit &, numerisch gleichen 
charakteristischen Zahlen jeweils zu linearen Elementarteilern gehoren, also (mehr- 
fache) Eigenwerte sind. €,... &, sind sodann Eigenwerte, x, ... %, sind die zugehi- 
rigen Eitgenvektoren. 


Die Treppeniteration konvergiert insbesondere (abgesehen von singuldren 
Fallen) fiir Hermitesches, in praxi fiir reell-symmetrisches A, soweit keine ent- 
gegengesetzt-gleichen Eigenwertpaare vorkommen. Falls A semidefinit ist, 
erledigt sich auch die letztere Einschrankung. Zu diskutieren ist noch die schon 
in 1.1 erwahnte Beschrankung fiir 7, genau formuliert die Bedingung / <1 —7, 
wenn A die charakteristische Zahl Null rv-fach besitzt. Offensichtlich lasst sich 
auch der noch fehlende Teilraum leicht ermitteln. Man kann nach der mit 
1=n—vy durchgefiihrten Iteration /’in beliebiger Weise zu einer vollen unteren 
Dreiecksmatrix /”’ erganzen, worauf sich ohne weiteres ein R’ bestimmen lasst 
derart, dass A J’ = J” R’ ist. 

Satz 1 ist ein Spezialfall des folgenden Satzes, aus dem hervorgeht, ob eine 
gewisse Spalte von J’; iiberhaupt konvergiert und wenn, mit welchem Konver- 
genzfaktor. 


Satz 2. Die k-te Spalte von I’; konvergiert (abgesehen von singuldren Fallen), 
falls 
lEn| > |€e41| oder falls §& ='.41:=/7 


und alle von y verschiedenen &, auch betragsmassig von y verschteden sind, und 
nur dann; das k-te Diagonalelement von R, strebt gegen &,. Im ersteren Fall 
herrscht lineare Konvergenz mit einem Konvergenzfaktor 


a Janel 
n= | E, 


? 


im letzteren Fall nur, falls alle charakteristischen Zahlen vom numerischen Wert y 
(mehrfache) Eigenwerte sind. Der Konvergenzfaktor ist dann 


by 


i a E, 


? 
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wobet &, die betragsgrisste unter den auf &y,44 folgenden (betragskleineren) charak- 
teristischen Zahlen 1st. 


v 
tigen J’, im ganzen es nicht tun. Darauf baut sich die im nachsten Abschnitt 


zu besprechende Methode zur Behandlung des Falles betragsgleicher charak- 
teristischer Zahlen (und héherer Elementarteiler) auf. Dieser Fall ist praktisch 
auch dann anzunehmen, wenn gewisse Spalten infolge eines nur wenig von Eins 
verschiedenen Konvergenzfaktors ungentigend konvergieren. 


Einzelne Spalten von J’; mégen also konvergieren, auch wenn die /-spal- 


2.3. Der Fall betragsgleicher charakteristischer Zahlen 


Dieser Fall stiitzt sich auf die Konvergenzaussage des folgenden Satzes 3. 
Mit I',(a, «) sei die aus den « Spalten a,a+1,...,a@+a-—1 von J’; gebildete 
rechteckige Matrix bezeichnet, mit T,(a, a) die aus den Zeilen a, a+ 1,..., 
a+a—1 von I',(a,«) gebildete « x a-reihige untere Dreiecksmatrix mit 
Diagonaleinsen. L,(a, «) sei definiert durch 


zw 


T(a, «) = I';_4(a, #) L,(a, a) . (10) 
Satz 312). Ist 
[Ee-al> [él wd [645-11 > lEsaol 2), 


so geht, abgesehen von singuldren Fallen, der von den o Spalten 
KiikeP Lo hd, RP Ole 


von I’; gebildete Teilraum im Limes i > 0o in sich diber; es gilt mit 


Talk o) —Pikpot baal 0 
sogar : 
Pylh, o) = TU VLA 0 eee 


Das bedeutet insbesondere, dass die / x J-reihige untere Dreiecksmatrix 
L,=L,(1, 1), die unter den Voraussetzungen von Satz 1 gegen die Einheits- 
matrix geht, unter den schwacheren Voraussetzungen von Satz 3 wenigstens 
unterhalb gewisser Ausbuchtungen folgender Art 


Ahnlich wie in Satz 2 sind auch in Satz 3 und den Zusdtzen offensichtlich schwadchere 
Voraussetzungen im Fall gleicher Eigenwerte méglich, auf deren explizite Formulierung an dieser 
Stelle wir glauben verzichten zu kénnen. 

13) Die erste bzw. zweite Ungleichung entfallt fiir k = 1 bzw. k to —1 =n. 


12) 
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im Limes 7 + co elementweise gegen Null geht. Dies findet, grob gesprochen, 
um so starker statt, je weiter das betreffende Element von der Diagonale 
entfernt ist. Es gilt naémlich 


Zusatz 1 zu Satz 3: Das Element in der j-ten Zeile und k-ten Spalte, 
k <7, von L,(1, l) strebt mit linearer Konvergenz gegen Null, falls 


1&.| > |& |, der Konvergenzfaktor ist a. 


ee Ad ve 


| 


Ein anschauliches Ergebnis liefert ferner der folgende 


Zusatz 2 zu Satz 3: Mit A,(1, 1’) wérde der l'-rethige linke obere Abschnitt 
von L, R; bezetchnet. Es strebt 


A I,(1, l’) — F,(1, 1’) Agai(1,l’) +0 (¢ +00), 


wenn Ll’ die Nummer irgendeiner honvergenten Spalte von I';, zum Beispiel der 
letzten konvergenten, ist). 


Es wird also dann A,(1, /’) im Limes 7 +o eine verkiirzte Ahnliche zu A, 
die wegen Zusatz 1 hochgradig reduzibel ist; jede nach Satz 2 konvergente 
Spalte erzeugt eine Einkerbung. Die charakteristischen Zahlen der (gering- 
dimensionalen) einzelnen Diagonalkastchen A,(k, o) streben gegen charakte- 
ristische Zahlen von A. Man beachte, dass die Folge der A,(1, 2’) und die 
der L,(1, 2’) dabei nicht zu konvergieren brauchen. Praktisches Interesse hat im 
allgemeinen nur der Fall einfacher Paare konjugiert-komplexer oder ent- 
gegengesetzt-gleicher Eigenwerte. Wir wollen ihn als Beispiel zur Erlauterung 
des vorangehenden ansehen: Als Abbrechkriterium fiir die Iteration diene nun 
etwa (mit den Bezeichnungen des vorigen Abschnitts) 


SN folic tar bs 100) 21 
) P(h, 1) —T cal’, 1) — ba Leak +1, 1)| Se 


fiir passend gewahlte /; , (J;,; = 0) und falls auf jedes nichtverschwindende /y, , 
ein praktisch verschwindendes /y ;,,, folgt. 
Es moge etwa ly, ,, nicht verschwinden. Ferner sei (ly 41, = / gesetzt) 


tL -0 


Lye 2)= (7 9) und Ryssbs2=(5 2) 


Dann hat 


a b 
Aha ime lag adie Riise Mele a yar 


die gesuchten charakteristischen Zahlen é,, &,.. 


14) Fiir / = nist l’ = 1 =n, es gilt dann trivialerweise sogar 
AI,(1,n) = £1, 2) Ay410, 2) - 
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Auch wenn es sich dabei um ein konjugiert-komplexes Paar handelt, kann 
man Real- und Imaginirteil der Eigenvektoren x,, x,,, in reeller Rechnung 
erhalten. Wir deuten einen Weg kurz an: Es sei 


‘ a ae A=a—ac, 
" toat Bp, wo a 
E cy B=V\A|Vsen4; 


fiir al +0 sind Eigenvektoren™) des Kastchens Ay ,,(%, 2) gegeben durch 


(2) =P (5). 


Sonach hat man einerseits mit 


y t= . 
Vy SB O= a Yuga 1; gga sree rea 


andererseits mit 
ye eat Me: vy = 0)8 Wo hg Sy = () 


“+1 n 


das Heraufrechnen an der Matrix Ay ,,(1,/) zu beginnen!®). Die Vektoren 
yn = {9, und y= {v4} bzw. die Vektoren z’= I'y,1 9’, ¥" = Ty 419" liefern 


x0 =x + Bx" Sy 


die Figenvektoren . von A zu den Eigenwerten : 
: tis = pat ow 


§ 3. Stabilisierte Simultaniteration 
3.1. Orthogonalisation an einem festen V ektorsystem 


Bereits einleitend erwahnten wir, dass die fiir die Treppeniteration kenn- 
zeichnende Dreieckszerlegung als Orthogonalisierungsprozess an den Zeilen der 
Einheitsmatrix E aufgefasst werden kann. Statt dessen kann auch eine allge- 
meinere nichtsingulare Matrix C benutzt werden, deren / erste Zeilen zur 
stabilisierenden Orthogonalisation dienen. 0; sei an C orthogonalisiert. Dann 
ergibt sich 0;,, durch Orthogonalisierung von A 6; an C, das heisst 


B= 6, Ry; (11) 
A 6; = 6544 Rigi, (12) 


wo R, wiederum eine obere Dreiecksmatrix und C 6, eine untere Dreiecks- 
are ist (die selbstverstandlich nicht explizit perce: zu werden braucht). 


15) Ein Eigenvektor, falls B = 0; dann ist i Hauptvektor. 


16) Wenn sich dabei der direkten Rechnung weitere Ausbuchtungen in den Weg stellen, ist 


allenfalls die triviale Lésung eines Systems von zwei Gleichungen mit zwei Unbekannten zu bewerk- 
stelligen. 
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Falls die Iteration konvergiert, also 0; > 0 strebt, kommt man sinngemass auf 
die Gleichung (3) bis (6) von § 1 zuriick. ist die an C orthogonalisierte Matrix 
(irgendwelcher) Eigen- (und Haupt-) Vektoren. Abgesehen davon, dass in der 
Rechenvorschrift A die Operation der Dreieckszerlegung durch die allgemeinere 
der Orthogonalisation an C [triang(B,) durch C-orth(B,)] zu ersetzen ist, 
andert sich nichts an der Bestimmung der Eigenwerte und Eigenvektoren. 

Es ist zu erwarten, dass auch die Konvergenzsatze im wesentlichen giiltig 
bleiben. Es gilt ja nach (11), (12) mit P, nach (8) . 


Att = 0.2; (wis) 
und, wenn 
CO. I", (14) 
gesetzt wird, die (7) genau entsprechende Dreieckszerlegung 
CatRe PP: (15) 
es ist dann 
OS AB Po! (16) 


In der Tat wird von der Wahl von C nur die Minorenbedingung (**) betroffen, 
und es ist jetzt immer méglich, durch geeignete Wahl von C (insbesondere als 
Permutationsmatrix, siehe unten) diese Bedingung ebenso zu erfillen, wie es 
durch die Wahl von B, fiir die Bedingung (*) bisher schon méglich war. 

Diese Satze 1, 2 und 3 von § 2 bleiben also, unter Abschwachung der For- 
derung beziiglich der auszuschliessenden singularen Falle, giiltig. 


3.2. Pivotelimination 


Insbesondere wird C Abschnitt einer Permutationsmatrix, wenn an B; 
anstatt der starren Dreieckszerlegung Pivotelimination durchgefiihrt wird, das 
heisst als Eliminationselement (pivot) das jeweils betragsgrésste unter den 
zulassigen Elementen (den Elementen, deren Zeilennummer noch keinen der 
vorhergehenden pivots gestellt hat) gewahlt wird. Die Pivotverteilung kann 
sich (und wird sich im allgemeinen) im Lauf der Iteration andern; da jedoch 
jedes Zwischenergebnis als Ausgangsnaherung angesehen werden kann, ist zur 
theoretischen Konvergenz nur zu fordern, dass schlzesslich eine feste Pivot- 
reihenfolge gewahlt wird. Vom numerischen Standpunkt ist dagegen nichts 
einzuwenden, es besteht kaum eine Gefahr, dass durch zu friihe Festhaltung 
der Sinn der Pivotelimination, ein unzulassiges Kleinwerden des Eliminations- 
elements zu verhindern, gestort wird. Im Gegenteil, man wird aus Aufwands- 
griinden bestrebt sein, méglichst frith die Pivotreihenfolge festzuhalten. 

Unter Ausserachtlassung der schliesslichen Festhaltung der Pivotreihen- 
folge besteht die Gefahr, dass es zu einer standigen Pendelung kommt. Wenig- 
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stens unmittelbar vor einer Konvergenzpriifung sollte also die Pivotreihenfolge 
festgehalten werden. | 


§ 4. Gegenseitig orthogonalisierte Bi-Iteration 


4.1. Rechenvorschrift der Bi-Iteration 


Als Nachteil durfte man bisher (wenigstens was die Eigenvektorbestimmung 
anbelangt) anfiihren, dass nicht die Matrix der Eigenvektoren direkt, sondern 
nur eine gewisse Zerlegung von ihr approximiert wird. Damit wird bei der 
nachtraglichen Berechnung der Eigenvektoren die erreichte Genauigkeit infolge 
von Rundungsfehlern in nicht ganz tibersichtlicher Weise gestort. 

Dieser Mangel erledigt sich merkwiirdigerweise von selbst bei der nachfol- 
gend diskutierten gegenseitig orthogonalisierten Bi-Iteration. 

Die beste Wahl fiir die Orthogonalisierungsmatrix C ware offensichtlich die 
Matrix der Linkseigenvektoren A. Beim Verfahren von Kocu [4] wird auch 
gerade durch Orthogonalisierung an der Matrix der schon bestimmten Links- 
eigenvektoren eine Stabilisierung in der Bestimmung der héheren Eigenwerte 
erreicht. Schon einleitend haben wir darauf hingewiesen, dass neben der bisher 
allein behandelten Rechtsiteration an Spaltenvektoren ebensogut eine Links- 
iteration an Zeilenvektoren moglich ist, die nunmehr (bei Transposition aller 
Formeln, jedoch nicht von A) Linkseigenvektoren liefern wiirde. Es legt 
nahe, beide Iterationsprozesse nebeneinander her laufen zu lassen, nunmehr 
jedoch die Iterationsspalten an der Matrix der Iterationszeilen und umgekehrt 
zu orthogonalisieren. Dabei bleibt das System der Iterationsspalten und das 
der Iterationszeilen standig gegenseitig orthogonal, und es ist zu erwarten, dass 
sie gegen ein System von Rechts- und von entsprechenden Linkseigenvektoren 
streben. Formelmassig haben wir, wenn die Matrix der / Iterationszeilen bzw. 
Spalten mit ®, bzw. Y; bezeichnet wird 


AW; =P 441 Rist, (17) 
@; A= Li 44 D544 , (18) 


wo R,;,, und L;,, obere bzw. untere Dreiecksmatrizen sind und die Zerlegung 
dadurch bestimmt ist, dass 


41% p= E (19) 
ist. Wir haben also folgenden Algorithmus: 


Rechenvorschrift B. Gegeben A und die Ausgangsmatrizen Cie Co Bie Bs 
von / Zeilen bzw. Spalten. 
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Fiir i = 0(1) N 
triang(C; B,) > ER, , 
=a 
B, R; SY, 


LC > By, 
A YY, a 0 Pe Be 
@; A ag OF 


¢ = N, falls 
|Pirr—ay P|] <<e und 


|®;,,— 6; ®;| <e fiir passende «,, Be; 


diag(®y ,, A Py 41) >{é,} ; 
Dy, >{m,}. 
Pryat > {to}: 


{é, }, {w, }, fas} o=1,...,/ sind beziehungsweise Eigenwerte, Links- und 
Rechtseigenvektoren. Die aus ®y., A Wy, bestimmten Eigenwerte haben 
Fehler von der Ordnung ¢ 7, wenn ®y ,, bzw. Wy ,, Fehler von der Ordnung ¢ 
bzw. 7 haben. Die Anzahl der genauen Stellen von &, ist also etwa die Summe 
der Anzahlen genauer Stellen von ®y ,, und Wy ,,. Fiir symmetrische Matrizen 
fallt natiirlich fiir C = B’ Links- und Rechtsiteration numerisch zusammen, 
hier verwendet man vorteilhaft Choleskiartige Zerlegung?’). Die Iterationsvek- 
toren ergeben sich dann von selbst auf Eins normiert. Sonst kann die Dreiecks- 
zerlegung beliebig normiert werden. Durch geeignete Normierung kann (und 
soll) die Gréssenordnung der einzelnen Zeilen von ®; und Spalten von ¥; stabil 
gehalten werden. 


4.2. Zusammenhang mit der Treppeniteration 


Analog zu (8) gilt 


At B= PRRs." “Ro (20) 
und 
CAt=L,£,---£,@,. (21) 
Daraus folgt wegen (19) 
CA2?#B=L,L,:---L,R,R,_,-:: Ry. (22) 
Vergleich mit § 3 ergibt, dass . 
| Oe Cae oe ey (23) 
und 
R; R;-1 ++: Ro = Pa: (24) 


17) Beachte, dass BY B, positiv semidefinit ist. 


ZAMP VIII/15 
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ist, wobei Pe und P,, diejenigen Matrizen sind, die sich [Gleichung (23)] bei 
stabilisierter Iteration mit der Orthogonalisierungsbasis C und der Ausgangs- 
naherung B ergeben. Wird noch in Verallgemeinerung von (10) 


L,=(C ra meal ee TS os (25) 

gesetzt und hat auch R,, die alte Bedeutung, so gilt 
LE, = Ly; -1 Li; (26) 
R; = Rs; Roj-1+ (27) 


Dies zeigt bereits, welchen engen und analytischen Zusammenhang die Bi- 
Iteration mit dem Friiheren hat. 
Aus 
GAP +l Bat Deady cael sds enol geen 


im Sinn der stabilisierten Iteration und 
CA2tA1 B= CANAALB = DPD, ++) 1gpOpAY RRs ae 


folgt 
OR Se Rogans (28) 
Ferner haben wir 


0, AY, = ®, Nr Rossa Reins = Loins Leics Pisa ar 


woraus folgt 
DM eT ts (29) 


2i+1 244+2 


und 


DD wee Pre 


a+1 a 2744+2 


(30) 


27+1-° 


Eine Reihe weiterer Beziehungen kénnte sich anschliessen, etwa nach (20) bzw. 
(21) sowie (22) 


CH Lyd honk, ie Rive sae Ry , (31) 
OD, Biz Ls Len <> Ligh eee (32) 


4.3. Konvergenzsdtze 


Fir die Konvergenzuntersuchung ist es zweckmassig, die ®; und WY, auf die 
Iterationsmatrizen 6; von § 3 zuriickzufiihren, nach Gleichung (16) definiert als 


0,, =A®? BPs. 


Aus (20) ergibt sich mittels (24) YW, = A‘ B P= = A- (A2* B P53), also die 
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grundlegende Beziehung 
Y= A-? Gs7 (33) 


In gleicher Weise lasst sich ®; auf eine an B stabilisierte Iteration an Zeilen- 
vektoren, mit der Ausgangsniherung C, zuriickfiihren. Aus Symmetriegriinden 
gentigt es, das Grenzverhalten von ¥Y; weiter zu diskutieren. In § 2 und § 3 
haben wir gesehen, dass 6,,, falls es konvergiert, gegen die an C orthogonali- 
sierte Matrix irgendwelcher Eigen- (und Haupt-) Vektoren {x} strebt, und 
dass im Fall der normalen Konvergenz die x, sich nach einer Ordnung der 
Betrage ihrer charakteristischen Zahlen einstellen. In der k-ten Spalte von 
6,; sind, sofern ihre Konvergenz linear ist, die Beitrage spiterer x,, 0 >, klein 
von der Ordnung 2 7 gegeniiber den Beitragen der ersten k x,. Nunmehr werden 
durch Anwendung von A~? die Beitrage derjenigen Eigenvektoren x,, die zu 
den betragskleinsten der ersten & charakteristischen Zahlen gehéren, derart 
angehoben, dass alle iibrigen x, ihnen gegeniiber klein von der Ordnung 7 
werden. Man kann daraus ersehen, wieso es der bistabilisierten Iteration etwa 
bei betragsverschiedenen charakteristischen Zahlen gelingt, die k-te Spalte 
von ¥; gegen die Richtung des k-ten Eigenvektors zu zwingen. 

Falls jedoch die k-te Spalte von 6,; nur logarithmisch konvergiert, schiesst 
die Anwendung von A~? weit iibers Ziel hinaus. Schon die einfachsten Beispiele 
zeigen, dass die k-te Spalte von ¥; nicht gegen einen Hauptvektor konvergieren 
kann?8). Wir werden den Fall héherer Elementarteiler (soweit er nicht, falls 
1 <n, fiir die Konvergenz der / Iterationsspalten bei normaler Konvergenz- 
tendenz belanglos bleibt), weiterhin nicht in Betracht ziehen. 

Damit diirfte verstandlich sein, dass die Konvergenzuntersuchungen, 
derentwegen wir auf den Anhang (A3) verweisen, auf gegentiber Satz 1 bis 3 
nur unwesentlich abgednderte Konvergenzsatze fiihren. Unter Konvergenz ist 
dabei Konvergenz der Richtungen der einzelnen Vektoren, das heisst Konver- 
genz bis auf Skalenfaktoren zu verstehen. Fiir eine Bedingung, unter der die 
Gréssenordnung der Vektoren stabil bleibt, vergleiche A3. In praxi wird es 
meist geniigen, die betragsgréssten Elemente der einzelnen Zeilen von Y; oder 
Spalten von @®; gelegentlich oder standig auf Eins zu normieren. 


Satz 4: Die k-te Zeile der D; und k-te Spalte der '¥,, Rechenvorschrift B, kon- 
vergiert (abgesehen von singuldéren Fallen) unter denselben Bedingungen, unter 
denen sie bei der Treppeniteration nach Satz 1 und 2 linear konvergiert. 

Die Konvergenz ist linear mit dem Konvergenzfaktor 


18) Die Divergenzerscheinung ist von der Art, dass zwar fiir jedes gentigend hohe t die Spalte 
einem (sich mit 7 4ndernden) Hauptvektor bh = b(t) beliebig nahe kommt, dass jedoch die Folge 
der h(i) gegen den zugehérigen Eigenvektor strebt. 
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wo &, bew. &, die &, betragsmdssig ndchsten der von betragsverschiedenen charak- 
teristischen Zahlen sind. Die Rechts- und Linkseigenvektoren treten geordnet nach 
den Betrigen ihrer Eigenwerte auf. Der Iterationsprozess ist numerisch stabil und 
selbsthorrigterend. 


Die Selbstkorrektur erhellt sofort daraus, dass die Iterationsvorschrift 
wieder eine Stationaritatsbedingung ist. Damit ist die Bi-Iteration eines der 
wenigen Verfahren, das «héhere» Eigenvektoren unmittelbar und selbstkorri- 
gierend liefert. Sie diirfte damit auch zur Nachiteration der Ergebnisse, die 
mit anderen, rascher konvergenten Verfahren zu erzielen sind, vorziiglich 
geeignet sein. 

Im iibrigen kénnen die einzelnen Spalten von VY; oder Zeilen von ®; wieder 
verschieden rasch konvergieren; wenn einzelne nicht konvergieren, braucht 
dadurch die Konvergenz der nachfolgenden nicht gestért zu werden. 

Der Fall betragsgleicher Eigenwerte erledigt sich diesmal wesentlich ein- 
facher. Die Matrix A,;,,= ©, A Y,=1,,;,, Ry,,1 zerfallt im Limes 7 > 0o in 
Kastchen (die die Eigenwerte gleichen Betrags enthalten), wenn die Iteration 
grossenordnungsstabil gefiihrt wird. 


Satz 51°): Unter den Bedingungen von Satz 3 wird bet groéssenordnungs- 
stabtler Iteration der von den Spalten k, k+1,...,k+o—1S1 von ¥, (ent- 
sprechendes gilt fiir ®,;) gebildete Teilraum durch A im Limes in sich transformiert, 


AW, (k, 0) — ©; (k, 0) Aga, (Rk, 0) >0 (i S00). 


Das Element in der j-ten Zeile und k-ten Spalte von Ag; ,, strebt mit linearer Kon- 
vergenz gegen Null, falls | €;| + | &, |, der Konvergenzfaktor ist 


vin (LB. 


4.4. Weitere Zusammenhange 


Falls /= n ist, gelten eine Reihe weiterer Beziehungen, die theoretische 
Zusammenhange mit der LR-Transformation von RUTISHAUSER [8], [9] und 
der AX-Transformation von RUTISHAUSER und dem Verfasser [7] aufzeigen. 

Da dann 0; L;.,= 6;,, ist, ergibt sich fiir die stabilisierte Iteration von § 3 

R,L,= Li41 Ryaa (34) 
die Grundbeziehung der LR-Transformation oder 
L; Ri = R;* Lia (35) 


t+1 


was man nunmehr mit den Formeln (29), (30) fiir die Bi-Iteration vergleiche. 


19) Anmerkung 1%) gilt sinngemass. 
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Wegen Ly * C = ®, = >" ergibt sich ferner aus (31), (32) 
Poa bee he a Ry} = 0, A, P*7? (36) 


i+1 t+2 


und entsprechend 
®, ms | oe : Sota as Er R; Ry of: 2 R, D, an ane P, ®, , (37) 
mit den Bezeichnungen, die in [7] verwendet wurden. Aus ®, Y%, = E ergibt sich 
FP Ay aial¥ PY (38) 


die Grundbeziehung der A2-Transformation. Gleichungen (36) und (37) lehren, 
wie man mittels der AX-Transformation unmittelbar die Eigenvektorbestim- 
mung durchfiihren kann [3]. Durch einfache Umformung mittels (34) oder auf 
Grund von (29), (30) erhalt man die Beziehungen 

Y; =P, L, Re" LI, Pi ae Ly i-4 Re; =P, i L, hee Datos Roya L 


24-1 


(39) 


24 


und 
O;— Ly; Ro jy Lyi-2 an Ly" R, ® = R,; Lyin Raja +++ Ry Ly" ®, (40) 


die ebenfalls in [3] zur direkten Eigenvektorbestimmung mit Hilfe der LR- 
Transformation vorgeschlagen wurden?°), 

Bei der Treppeniteration und bei der Bi-Iteration ist man jedoch nicht an 
den Fall / =m gebunden. Bei Matrizen hohen Grades, wenn einige wenige 
Eigenwerte und eventuell Eigenvektoren zu bestimmen sind, diirften sich 
daher diese Verfahren empfehlen. 


4.5. Mischtypen von Bi-Iteration und LR-Transformation 


Offensichtlich kann jederzeit von der nicht-selbstkorrigierenden L R-Trans- 
formation mit Eigenvektorapproximation nach (39), (40) auf die selbstkorri- 
gierende Bi-Iteration tibergegangen werden und umgekehrt. 

Dadurch ist unter Umstanden (bei besonderer Gestalt der L,, R,) eine 
erhebliche Einsparung an Rechenaufwand erzielbar, ohne dass die Selbstkor- 
rektur prinzipiell aufgegeben werden muss”). Wir glauben, dass in praxi auch 
ein Mischtyp solcher Art die optimale Lésung hinsichtlich Sicherheit und 
Aufwand darstellen kann. : 


20) Fiir den Spezialfall, dass By = (t, At, ..., A-1t) und Cg = (5, 45,...,A%15) ist und 
damit in Ly und Ry nur die erste untere bzw. obere Nebendiagonale nicht verschwindet, ergibt sich, 
falls 1 =n, die Rechenvorschrift B des QD-Algorithmus von RuTisHAUSER [10]. 

21) Falls J < n, sind jedoch jeweils nur die ersten 7 — B Zeilen von ®,, , bzw. Spalten von ¥, 
mittels der Formeln (39), (40) berechenbar, wenn in Ly und Ry, nur die ersten B unteren bzw. oberen 
Nebendiagonalen nicht verschwinden. Dies ist etwa der Fall, wenn A von Bandgestalt ist und 
By = Cy = E. Wie beim Spezialfall des QD-Algorithmus [10], Rechenvorschrift D (6B = 1), mtissen 
die B letzten Zeilen bzw. Spalten durch Iteration ermittelt werden. 
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Anhang 


A 1. Determinantendarstellung der Dreieckszerlegung 


Die Determinantendarstellung einer Dreieckszerlegung G = U DO, wo G 
cine nichtsingulire quadratische, U bzw. O eine untere bzw. obere Dreiecks- 
matrix und D eine Diagonalmatrix ist, lautet bekanntlich 


1 5 07eR aie 
(UO) jn = Vin oh Ue hits gto (0) in = Vin ee 
wo 
Su S12 Sik 
§21 S22 Sar 
Vin = ViplG@) = | eevee esc ne eee acne en aes F (41) 
Seti beatae so Spee | 
Sj1 §j2 -++ SiR 


Yo=l, Vik = Ypj(G*) 


und 0;, 04 beliebige von Null verschiedene Normierungsfaktoren sind. Ge- 
brauchlich sind die folgenden Normierungsvorschriften, wobei D = E wird: 


a) O,=Yews O% ~ Vi=i, faa? 
b) 0, = VBA bi 0; eae 


wodurch U bzw. O Diagonaleinsen bekommen, sowie fiir positiv-definites, 
hermitesches G , 


c) 0, =V7i. aes Reales 0; =Vyii Maia , 


wodurch U? = O wird. 

Falls G eine rechteckige Matrix von m Zeilen und / Spalten ist, die den 
Héchstrang hat, gelten offensichtlich die obigen Formeln (bei sinngemassem 
Indexlauf) auch fiir Zerlegungen, bei denen 


(i) U (fiir m > 1) Trapezgestalt bzw. 
(ui) O (fiir m <1) Trapezgestalt hat. 


A 2. Konvergenz bei Treppeniteration und der stabilisierten Iteration 


Nur zu Beweiszwecken sei die Matrix C von § 3 durch Hinzufiigen weiterer 
linear unabhangiger Zeilen zu einer quadratischen (nichtsinguléren) Matrix 
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C erganzt gedacht. Die Rechenvorschrift wird dadurch nicht berithrt, und wir 


werden darauf zu achten haben, dass auch itn Endresultat nur die Zeilen von 
C auftreten. 


Wir haben sodann die Dreieckszerlegung von CA‘ Bau diskutieren, 
G=CA‘B=T,P,, 


wobei C = E fiir die Treppeniteration. Dabei hat J’; Diagonaleinsen, wir haben 
also den Fall a) von Abschnitt 1 vor uns mit 


Wh 
(Li)j2 = ott 
kk 
Falls 1 < n, haben wir ausserdem den Fall (7), ’; von Trapezgestalt. Als Ele- 


mente von J’; (und P;) treten damit Minoren von C A‘B auf, deren Grenzver- 
halten schon PERRON [6] und Mintz [5] studiert haben. Fiir unsere Zwecke 
ergibt sich folgendes?¥) : 

A habe die kanonische Darstellung 


A= EW (WAX = BY. 
Dann haben wir 
Gat, 
mit den Matrizen U = CX, V = WB, gebildet aus den Gewichten fiir die Ent- 


wicklung der C, B nach dem System der Zeilen bzw. Spalten von W bzw. X. 
Wie-in einer friiheren Arbeit [2] ergibt sich damit fiir die k-te Spalte 
(k SJ) (L), von 1; 


Vik 
Ney Une. Us cp 
U2, U2 4, woe Ue 
yp = aaes Eat Cee een ie Pale ane, Lahm (42) 
1B be | Meta, “e—t,ay Uy, cx By B2--* BR 1 k 
| Uy, Uy, te, 


wo u,, die -te Spalte von U und 


LN py ag ay 


¥y VQ°8* MR 


der Minor aus den Zeilen py fa +** “4, und den Spalten », % --- ¥, von FF ist. 


22) Wegen des engen theoretischen Zuasammenhangs der Treppeniteration mit der LR-Trans- 
formation von RuTIsHAUSER geht unsere Darstellung in diesem Abschnitt weitgehend einem 
Beweisgang von RUTISHAUSER [9] parallel. 
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Nun ist 6, = C-1 J’, also ergibt sich wegen X = C-! U fiir die k-te Spalte 


(0;), ae (i) oY , 
Vik 
14, Uy x5 Uy ap | 
Ue op U2, vee Une | 
elas CSN VG. Gs Boa 5 Cae ee | (Anes, (lee, 
eeeth |e Mees i, es ee a 
Fay Fa Xap 


x,, ist die u-te Spalte von X, der Matrix der Eigen- (und Haupt-) Vektoren von 
A. Ferner ist in (42) und (43) explizit 


Vek = oe LU] 12. +k Ll eect, babpeare eee 


sae 4s Oy hg? Op By Bo-+ BR 1 aiag 
By Bg:-- 


Man beachte, dass in Gleichung (43) fiir (6,),, tatsachlich nur noch die ersten 


k <1 Zeilen von C eingehen, also diejenigen Zeilen von C, die in die Iteration 
eingreifen. 

Die Diagonalelemente von K — die charakteristischen Zahlen — seien nun 
nach der Ungleichung (9) geordnet. Von normaler Konvergenztendenz sprechen 
wir, wenn unter den Gewichten [Vlosps... Bt derjenigen Terme 


BI ee oe 
By Bo-++ Br 


die starker anwachsen als alle tibrigen Terme, wenigstens eines nicht ver- 
schwindet. 

Falls €,...&€ und alle mit é, numerisch gleich charakteristischen Zahlen 
Eigenwerte sind, ist dafiir hinreichend, dass 


LVivevmcy AOR (0 aaa Oe ced ee (*) 


1Q<.:k 


Die explizite Formulierung der etwas weiteren notwendigen Bedingungen fiir 
normale Konvergenz, insbesondere im Falle nichtlinearer Elementarteiler?’), 
interessiert kaum angesichts des Umstandes, dass es sich dabei offensichtlich 
um Bedingungen handelt, die gleich (*) nur in singuldren Fallen der Wahl von 
B verletzt sind, wenn man die Randbedingung fiir / (vgl. den letzten Absatz 
von 2.2) beachtet. 


23) Siehe zum Beispiel [2]. 
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Zur Konvergenz ist sodann noch zu verlangen, dass y\) nicht schwacher 


wachst als die tibrigen Elemente von ( bzw. 00), Unter den oben formulierten 
Bedingungen ist dafiir hinreichend, dass auch 


AP oe OS peo oe eA (=) 
12-.+k 

Auch hier tritt, ahnlich wie oben, nur in singularen Fallen der Wahl von C eine 
Konvergenzstérung auf. 

Die Aussagen der Satze 1 bis 3 sind nunmehr aus der Darstellung (43) 
(wobei im Fall der Treppeniteration U = X zu setzen ist) klar zu entnehmen. 
Fir Satz 3 ist lediglich zu bedenken, dass im Limes 7 > oo offensichtlich (6;) ,’ 
bzw. (1’;)v, k Sk’ Sk+c, in einem o-dimensionalen Raum bleibt. 


A 3. Konvergenz der Bi-Iteration 
Es wird geniigen, die Konvergenz von 
P,=A-*6,; 


zu diskutieren, wo 6, wie in Abschnitt 2 definiert ist. Ergebnisse fiir ®; ergeben 
sich dann durch Transponierung und Vertauschung von B und C. 

Aus (43) gewinnen wir sofort fiir die k-te Spalte (W,), von YW; (wobei wir 
y?9 durch 0( ersetzen, entsprechend der nunmehr offen gelassenen Diagonal- 


normierung) : grime, 
; (Ps) iat AG ne , 
a 
k 
bs fy U1 2 Wy x, 
Ue, Ue a ss Weg, t 
ayy | Qi 
oe | AT oe eR esl Gla ber LA Mlrar ee eae od 
Oy Ay + Op * By Bg: BR 12-..% 
By Ba++ Be Un 1,04 ULI a5, k—1, up 
= ie =i 
Al teat Xe At 2S ig 


Setzen wir den Fall linearer Elementarteiler voraus, so vereinfacht sich diese 
Beziehung zu 


M1 a4 1 a5 “1 ap i 
| 43 oy Yo09 42 op 
ee eee |! Sivetat He -aat le) (o's Seto el 0 nie oi tet ol tye 0) fae , 2% 
9) Ah LE jails inhib SS) clea agate’ 
u LBp><f 
Gy cay | Up aig, Ue-trag °° “h—lap 
1 (45) 

it i 3 
Xa, hy a2 e aR | 
ay ap Lp 
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Bedingung (*) verbleibt als Bedingung fiir normale Konvergenztendenz von 
Y;, Bedingung (**) kommt hinzu als Bedingung fiir normale Konvergenz- 
tendenz von ®;. Ist nun |&,_4|> |&,| = |v4i| = -** |Sn40-a|  |Sxte 
wachsen offensichtlich fiirx < k < «+ o— 1die zu den Vektoren x, bis x, ,,_; 
gehorigen Faktoren starker als alle itbrigen, die Richtung von YW; bleibt also im 
Limes 7 -+ oo in dem von x,.,,_1°*: %, gebildeten o-dimensionalen Raum. 

Insoweit sind auch die Konvergenzaussagen der Satze 4 und 5 offensicht- 
lich. Wenn pauper der k-ten Spalte von WY; herrscht, geht 
offensichtlich 0 wie |& ... &|" |& ...&_1|". Das bedeutet, dass das’ k-te 
Diswotateieent von J',; mit |é,|’ geht (anstatt Eins zu sein, wie bei der 
Treppeniteration). Mit diesem Resultat ergeben sich die Konvergenzfaktoren 
von Satz 5 unmittelbar. Eine andere Konsequenz der Gréssenordnungsstabili- 
tat ist, dass das k-te Diagonalelement von L; (und von R,) gegen |&,| geht. 
Man erreicht eine gréssenordnungsstabile Zerlegung offensichtlich, wenn die 
Dreieckszerlegung in der Rechenvorschrift B betragsmassig der Vorschrift c) 
von Abschnitt 1 entspricht, insbesondere also durch Choleskyartige Zerlegung 
im Falle hermitescher A. 

Im iibrigen andert sich, wie man leicht sieht, an diesen Aussagen auch im 
Falle héherer Elementarteiler nichts, wenn wenigstens &, bis €, und alle mit &, 
betragsgleichen charakteristischen Zahlen Eigenwerte sind. 


> 
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Summary 


We have described two closely related numerical methods for the solution 
of eigenvalue problems called ‘stabilized iteration’ and ‘ mutually stabilized bi- 
iteration’. The iteration process shows numerical stability, has no loss of leading 
figures and is self-correcting, giving some eigenvectors immediately in the second 
case, in triangular decomposition in an important special case of the stabilized 
iteration, called: ‘Treppeniteration’, 

The processes are mathematically closely related to RUTISHAUSER’s LR- 
Transformation, but differ in numerical aspect. A certain mixed process, based 
on RUTISHAUSER and the author's ideas, seems to be the optimum way with 
regard to both security and numerical labor. 


(Eingegangen: 27. November 1956.) 


Le décrochage tournant - 
des machines axiales génératrices monoétages. 
Lois de similitude 
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1. Introduction 
1.1. Nature du décrochage tournant 


La naissance du décrochage tournant est provoquée par le décollement de 
la couche limite sur l’extrados d’un ou de plusieurs aubages. 

On sait que, l’incidence du fonctionnement d’un profil augmentant, l’énergie 
cinétique de la couche limite finit par n’étre plus suffisante pour vaincre le 
gradient de pression positif de l’extrados: la couche limite décroche, ce qui se 
traduit par une chute brusque de la portance (figure 1). 


Direction de prapagation 
du décrochage tournant 


Figure 1 Figure 2 
Courbe du coefficient de portance en fonction Représentation schématique du décollement 
de langle d’incidence pour un profil isolé. dans une grille d’aubes. 


L’incidence de décrochage dépend de divers facteurs (nombre de REYNOLDs, 
turbulence...). Aprés décrochage, en diminuant l’incidence, on constate que 
le phénoméne a une hystérése [1]?). 

Considérons une grille d’ailes décélératrice (agissant en compresseur): 
lorsqu’on atteint l’incidence limite, une ou plusieurs ailes contigues décrochent ; 
une quelconque dissymétrie dans l’écoulement ou dans la construction de la 
grille suffit a provoquer ce décrochage préférentiel. La zone décrochée offre 
au flux une résistance augmentée: son débit est réduit ou méme annulé (blo- ’ 
quage du débit). Le courant incident est dévié de part et d’autre de la zone de 
décollement. Si les aubes voisines sont suffisamment proches, l’angle d’inci- 
dence est diminué du cété des aubes 1 et 2 et augmenté du cédté de l’aube 4 
(figure 2). Pour celle-ci, l’angle d’attaque dépasse la limite et le décrochage se 
produit, tandis que du cété de l’aube 2, la diminution de l’incidence provoque 
le réattachement de la couche limite. Ainsi dans une grille d’aubes d’extension 


*) Les chiffres entre crochets renvoient a la Bibliographie, page 249. 
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infinie, la zone décrochée se propage le long du front de la grille. La vitesse de 
propagation de ce phénoméne fait l’objet d’investigations théoriques et expé- 
rimentales actuelles [2], [3]. 

Considérons a présent un étage de machine axiale génératrice (ventilateur, 
soutflante, compresseur axial). A chaque rayon, rotor et stator constituent 
une grille d’aubes. Le phénoméne qui vient d’étre décrit pour une grille s’y 
produit donc lorsqu’on augmente suffisamment l’angle d’incidence des profils, 
c’est-a-dire lorsqu’on diminue le débit de l’étage considéré. 

Suivant la hauteur relative des aubes et la conception aérodynamique de 
Pailetage, l’incidence limite peut étre atteinte localement, par exemple au 
sommet ou au pied des aubages, ou elle peut étre réalisée quasi simultanément 
sur toute la hauteur des aubes. Suivant le cas, le décrochage se manifestera sur 
une partie seulement de la hauteur des ailettes ou sur toute leur hauteur; il 
peut apparaitre au rotor ou au stator. 

Deux types de décrochage tournant doivent étre distingués: le décrochage 
progressif et le décrochage brusque de la base au sommet. 


1.2. Le décrochage progressif 


Les propriétés principales du décrochage progressif sont les suivantes: 

La courbe caractéristique de l’étage (coefficient de pression y en fonction 
du coefficient de débit ®) est continue (figure 3). Aussi le décrochage est-il 
qualifié de progressif. Dans la zone de décrochage, y diminue avec ®. 

Il existe un et souvent plusieurs secteurs décrochés, disposés symétrique- 
ment. En général ils n’occupent qu’une partie de la hauteur des ailettes, c’est 
pourquoi le décrochage progressif est souvent qualifié en outre de partiel*). 


Sa 


$ 


Figure 3 


Courbe caractéristique: coefficient de pression y en fonction du coefficient de débit ®, pour un 
rotor avec décrochage progressif. 


3) Le fait d’étre partiel n’est pas le caractére principal du décrochage progressif. Le décrochage 
progressif peut en effet s’étendre a toute la hauteur des aubes. Il continue a étre, en général, formé 
de plusieurs zones et 4 avoir une répartition progressive des fluctuations de la vitesse axiale [5]. 
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Le décrochage progressif et partiel se produit donc sur les rotors a rapport de 
moyeu4) faible, pour lesquels une diminution du débit augmente plus rapide- 
ment les incidences en une région des aubes. Cette région peut étre le sommet, 
le pied ou méme une région intermédiaire, suivant la conception aé¢rodyna- 
mique de l’aubage considéré. 

L’extension tangentielle des zones décrochées varie avec le rayon, de 
méme que la répartition des fluctuations de la vitesse axiale AU/U (figure 4). 


AU 
U 


Zone décrochée temps 


Figure 4 
Représentation schématique du décrochage Fluctuations relatives de vitesse en fonction 
progressif lorsque celui-ci se produit ausommet . du rayon: 1 sommet des aubages; 2 dans la 
des aubes (les zones de décollement sont zone décrochée; 3 pied de la zone décrochée; 
ombrées). 4 pied de l’aubage. 


Les zones décrochées tournent avec une vitesse angulaire propre qui est 
une de leurs caractéristiques. Pour un rotor, elles tournent dans le sens de 
rotation de la machine avec une vitesse angulaire absolue proportionnelle mais 
inférieure a la vitesse du rotor. 

Une diminution du débit provoque une extension radiale et tangentielle 
des secteurs décrochés [4], [5], éventuellement un changement du nombre de 
zones de décrochage s’accompagnant d’une variation de leur vitesse angulaire 
((2], [3], [4], [5], [6], [7]). On a marqué a titre indicatif sur la caractéristique 
y — ® (figure 3) le nombre de secteurs décrochés. 


1.3. Décrochage brusque de la base au sommet 


Les propriétés de ce type de décrochage sont les suivantes: 

La courbe caractéristique de |’étage y — ® comporte une discontinuité a 
l’apparition du décrochage (figure 5); c’est pour cette raison que le décrochage 
est qualifié de brusque. 

Il existe un et parfois plusieurs secteurs décrochés symétriques [8] s’étendant 
toujours de la base au sommet des aubes. Le décrochage brusque se produit 
lorsque Vincidence limite est atteinte quasi simultanément sur toute la hauteur 
des ailettes. Cette circonstance est courante lorsque le rapport de moyeu est élevé. 


4) Rapport du diamétre du moyeu au diamétre extérieur de l’ailetage. 
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Pour un rotor, la vitesse absolue de rotation du décrochage est inférieure 
a celle de la machine. Comme nous le montreronis (2.2), elle ne dépend guére 
du point de fonctionnement dans la zone décrochée. . 

Le phénoméne présente une hystérése: le décrochage apparait A débit 
descendant pour une valeur du débit plus faible que celle le faisant disparaitre 
a débit croissant (figure 5). Cette particularité est causée par l’hystérése de la 
courbe de la portance en fonction de l’incidence (figure 1). Le décrochage 


> 


Figure 5 
Courbe caractéristique: coefficient de pression y en fonction du coefficient de débit ®, pour un 
rotor avec décrochage brusque. 


brusque peut étre précédé d’un décrochage progressif s’étendant sur une 
portion de la hauteur des aubes ou sur la totalité de cette hauteur [5]. 


1.4. Distinction avec le pompage 


Il convient de ne pas confondre le décrochage tournant avec le pompage. 

Le décrochage tournant est un phénomeéne qui affecte localement l’écou- 
lement. L’écoulement absolu est pulsé, mats V écoulement relatif (par rapport a 
des axes tournant avec la vitesse angulaire du décrochage) est un écoulement 
asymétrique permanent. La fréquence des pulsations est proportionnelle a la 
vitesse de rotation de la machine. 

L’effet du décrochage tournant reste localisé au voisinage du rotor et 
n’affecte pas le circuit. Les mesures de débit et de pression suffisamment loin 
en amont et en aval de la roue, mesures nécessaires pour déterminer les coeffi- 
cients y et ®, s’effectuent aisément, en écoulement non pulsé. 

Le pompage au contraire affecte l’ensemble de l’écoulement et le rend vio- 
lemment pulsatoire. L’écoulement absolu comme 1’écoulement relatif ne sont 
plus permanents. Le pompage affecte l’écoulement dans tout le circuit. La 
fréquence du phénoméne ne dépend pas de la vitesse de rotation de la machine. 

L’appellation pompage tournant ou rotatif parfois donné au décrochage 
tournant nous semble ambigué. 
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2. Propriétés de similitude 


Du fait des caractéristiques des rotors étudiés, nos investigations expéri- 
mentales ont porté uniquement sur le décrochage brusque du pied a la téte des aubes. 
Nos conclusions se rapportent uniquement a ce type de décrochage. 

Nos essais ont porté sur deux rotors caractérisés par des profils vrillés, un 
rapport de moyeu de 0,63 pour l’un et 0,76 pour l’autre. Le nombre d’aubages 
est respectivement de 14 et de 9. Les essais rapportés dans ce travail sont 
relatifs au rotor dont le rapport de moyeu est 0,63 et dont le rapport profon- 
deur/pas varie de 0,41 au sommet a 1,16 a la base. 


2.1. Similitude pour. le fonctionnement avec décrochage tournant brusque 


On trouvera a la figure 6 la courbe caractéristique relevée en fonctionnement 
normal et dans la région de décrochage, pour plusieurs vitesses de rotation. 
La courbe a été portée en coefficients sans dimensions: coefficient de pression 
y en fonction du coefficient de débit ® avec 
Zo ed U 


a ab ste Om 
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Figure 6 — 


Courbe caractéristique du rotor: coefficient de pression py en fonction du coefficient de débit ® 
relevés pour diverses vitesses, A débits croissant et décroissant. 
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fluide; 7,, le rayon moyen du rotor défini par 


ai ae 
oes |/ Yext -- Vint 
mn ee 


o) 


4 


On constate que les lois de similitude classiques sont aussi bien vérifiées dans 
la région avec décrochage qu’en fonctionnement normal, dans les limites ou 
Vinfluence de la variation du nombre de REYNOLDS avec la vitesse est négligeable. 

Le diagramme fait apparaitre l’hystérése du phénoméne de décrochage 
tournant. Lorsqu’on opére trés progressivement (pour éviter les effets dyna- 
miques) on observe que l’apparition et la disparition du décrochage s’ effectue 
toujours pour la méme ouverture de fonctionnement, quelle que soit la vitesse. 

Pour le fonctionnement avec décrochage, les lois de similitude sont évi- 
demment vérifiées pour les secteurs non décrochés. Nous venons de voir 
qu’elles sont vérifiées pour l’écoulement moyen. On peut en conclure qu’elles 
gouvernent également les caractéristiques propres aux secteurs décrochés. 


2.2. Lot de similitude dela vitesse du décrochage tournant brusque 


La vitesse de rotation absolue du décrochage tournant a été déterminée: 

a) pour diverses ouvertures de fonctionnement ; 

b) pour diverses vitesses de rotation du rotor; 

c) par différentes méthodes (tachymétre a lame vibrante, enregistrement de la 
pression totale instantanée devant et derriére le rotor, analyse spectrale du 
bruit et des vibrations). 

On trouvera a la figure 7 le rapport Np/N de la vitesse du décrochage 

tournant a la vitesse de rotation du rotor, en fonction du coefficient de débit B. 


es points expérimentaux a débit progressif croissant 
—o—e— points expérimentaux a débit progressif décroissant 


2|z 
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Figure 7 


Vitesse réduite du décrochage tournant déterminée par lame vibrante en fonction du coefficient 
de débit (points expérimentaux a débit croissant et décroissant). 


ZAMP VIII/16 


242 Stmon GOLDSTEIN et ANDRE L. JAUMOTTE , ZAMP 


On constate que pour un coefficient @ fixé, c’est-a-dire pour une ouverture 
donnée, le rapport Np/N est constant. Pour une ouverture de fonctionnement 
déterminée, la vitesse de rotation du décrochage tournant est proportionnelle a la 
vitesse de rotation de la machine. Cet énoncé constitue la loi de similitude du 
décrochage tournant brusque. 

La méme figure 7 fait apparaitre une propriété du décrochage tournant 
brusque, signalée en 1.3. 

La vitesse absolue du décrochage tournant ne dépend guére du point de 
fonctionnement. Elle se reléve légérement prés de la disparition du décrochage 
et au voisinage de l’apparition de l’écoulement avec courant de retour (faibles 
ouvertures). 

Pour le rotor essayé, on obtient 0,48 comme valeur du rapport Np/N. 


3. Détermination expérimentale de la vitesse de rotation 
du décrochage tournant 


Le décrochage tournant est caractérisé, pour des conditions d’amont uni- 
formes, par un écoulement absolu asymétrique dans l’espace et pulsé dans le 
temps, au voisinage du rotor. 

Les fluctuations au voisinage du rotor (en amont comme en aval) affectent 
localement la vitesse du fluide, sa pression et sa température. Elles induisent 
des vibrations et un effet sonore. 

On pourra donc en principe déterminer la vitesse de rotation du décrochage 
tournant par tout dispositif sensible aux caractéristiques locales instantanées 
de l’écoulement (pression, vitesse, température) ou aux effets induits (bruits, 
vibrations). 


3.1. Détermination par une mesure dans l écoulement 


3.11. Lame vibrante et stroboscopie de fils. Un tachymétre a4 lame vibrante de 
longueur réglable est placé dans l’écoulement avec la lame dirigée sensiblement 
dans la direction générale de l’écoulement et placée de fagon a vibrer dans la 
direction périphérique. 

La longueur de la lame est réglée pour obtenir la résonance avec le phéno- 
méne étudié, c’est-a-dire le décrochage tournant. La fréquence de celui-ci est 
lue directement a l'appareil. La résonance est nettement marqueée. © 

Notons que les perturbations dues au passage des ailettes sont trop faibles 
pour permettre de déceler une résonance a une fréquence égale au produit de la 
fréquence de rotation de la roue et du nombre d’aubes. 

Nos mesures, effectuées dans la gamme des fréquences de 10 4 20 Hz nous 
ont montré qu’avec le tachymétre utilisé, la + urecinieg était de l’ordre de +2%. 
La mesure est simple et rapide. 


Par déplacement radial de l’appareil, il est possible de déceler si le décrochage 
s’étend sur toute la hauteur de l’ailetage. 
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La stroboscopie de fils (par exemple une rangée de fils sur un support 
placé radialement) est basée sur le fait qu’un fil fin suit les fluctuations de 
directions de la vitesse. La fréquence est détérminée par stroboscopie. 

La méthode est plus difficile 4 mettre en ceuvre car, contrairement A la 


lame vibrante, le fil est sensible aux fluctuations radiales comme aux pulsations 
tangentielles. 


La précision est faible. 

Notons que les deux méthodes signalées ne permettent pas de déceler le 
nombre de zones décrochées. En fait, elles donnent le produit de la fréquence 
de rotation par le nombre de zones décrochées. 


3.12. Pression, vitesse ou température instantanées. Une sonde de pression 
statique ou de pression d’arrét placée en amont ou en aval du rotor et raccordée 
a un capteur de pression dynamique permet la détermination de la fréquence 
de rotation du décrochage rotatif. 

Nous avons opéré avec une sonde de pression d’arrét en amont et en aval 
du rotor a divers rayons, avec une prise de pression statique, en amont et en 
aval a la paroi. 

Le capteur de pression utilisé est du type capacitif. La fréquence peut étre 
déterminée par analyse harmonique de la tension du capteur ou par enregis- 
trement de cette tension. 

On trouvera a la figure 8 des enregistrements de la pression d’arrét au rayon 
relatif 0,85, immédiatement en amont et en aval du rotor pour deux vitesses 
de rotation. Ces divers enregistrements ont été effectués avec la méme ampli- 
fication: ils sont done directement comparables. Ils ont été effectués avec un 
enregistreur 4 plume répondant jusqu’a 100 c./s. On constatera que les fluctua- 
tions d’aval sont les plus importantes. 

D’une maniére tout-a-fait analogue, les fluctuations de vitesse peuvent étre 
décelées par anémométre a fil chaud. Nous n’avons pas utilisé cette technique 
a l'Institut de mécanique appliquée de l'Université. 

Au lieu d’une sonde de pression ou de vitesse, on peut utiliser une sonde de 
température. Cette méthode a été appliquée avec succés par E. L. CostiLow 
et M. C. Huppert [5]. 

Les mesures par pression, vitesse ou température instantanées permettent 
de déterminer la vitesse de rotation du décrochage et le nombre de zones pourvu 
que l’on enregistre simultanément les indications de deux capteurs décalés l'un 
par rapport a l’autre d’un angle que l’on peut faire varier. 


3.2. Détermination par un effet indut 


A notre connaissance, la détermination de la vitesse de rotation du décro- 
chage tournant par un phénoméne induit comme le bruit et les vibrations n a 
jamais été employée. Elle n’est d’ailleurs utilisable que pour un rotor isolé. 
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Investigation au rayon 7/7, = 0,85 


; + Base de temps 
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Figure 8 


Enregistrements filtrés de la pression d’arrét au rayon relatif 0,85, en amont et en aval du rotor, 
a deux ouvertures différentes, et pour deux vitesses de rotation. 


245 


trices 


‘ 


énéra 


Décrochage tournant des machines axiales g 


Vol. VIII, 1957 


‘uol1e101 ap souenbary “f fasvyoos0ap op souenbgay Pf *(u/"1 QELT :9}UeJSUOD IOjOI Np assez) 
jueuIMo} aseyooro9p np uorpuedde seide yo yueAe 10}OI Np jMIq Np jefduroo enbruourrey ssAjeuy 


6 oInsly 
ZH) OOv OO€ 002 ool 
Seley ] a i a, i, al a ae a al 
| | | 
: |; fp nf ak ee 
Ae | aS a 
| u eS ) | | 
Tee sage 7a Tie tke ee is 
> | | 4 “fh | i] ws | ir 4 
! | 
a me = Pa. \ Ty [ | | | i 
eet e gt Wd Te ce SAL | La 
Scare <i i. ae j ; 
\ . || | [ 
poe ee a + po} 4 ++ + on ; 
\ i 
\ | ; é 
7 +} jfedat H - + 9é¢'0 ig } afeydasap saide advanbasj ap aiyzadsg---+- ; ad 
\ / i oveo= 
2 - 
A Se + a | IseO= cb aBeusoisap JueAe aduvanbaJ) ap asjgadg = }—4—— 
) aE TaN zled= > 


246 Simon GOLDSTEIN et ANDRE L. JAUMOTTE ZAMP 


3.21. Analyse spectrale du bruit de la machine. L’apparition du décrochage 
tournant est marquée par une modification du bruit émis par la machine. 

La méthode consiste 4 comparer le spectre de fréquence du bruit avant dé- 
crochage au spectre du bruit mesuré pour un fonctionnement avec décrochage. 

La comparaison des deux spectrogrammes révéle immédiatement le pic 
d’intensité apparaissant avec le décrochage tournant: la fréquence corres- 
pondant au sommet de ce pic est la fréquence du décrochage tournant ou 
plutdt le produit de cette fréquence par le nombre de zones décrochées. 

La méthode est apparue trés précise (+ 0,5% entre 10 et 20 Hz). Elle ne 
permet pas de déterminer le nombre de zones de décrochage. 

Il faut noter que le bruit n’est pas d’intensité constante; la détermination 
de l’intensité moyenne en fonction de la fréquence, a l’analyseur harmonique, 
exige pour chaque point une moyenne de 10 a 25 lectures effectuées a des 
intervalles de temps égaux (par exemple 3 s). 

On trouvera a la figure 9 analyse harmonique complete du bruit d’un 
ventilateur axial, avant et aprés décrochage, jusqu’a une fréquence de 500 Hz. 

La figure 10 donne le détail de cette analyse, au voisinage de la fréquence 
du décrochage tournant, effectuée avant et aprés décrochage, a une ouverture 
déterminée et a diverses vitesses de rotation. Le rapport K est celui de la 
vitesse absolue du décrochage tournant a la vitesse de rotation de la machine. 


3.22. Analyse spectrale des vibrations. Un capteur de vibrations est fixé sur 
l’enveloppe du ventilateur (figure 11) ou sur son moteur d’entrainement, s’il 
n’y a pas d’accouplement élastique. 

Le capteur de vibrations utilisé donne une tension proportionnelle a la 
vitesse de déplacement du corps vibrant. 

L’analyse de cette tension (a l’analyseur harmonique) pour un fonctionne- 
ment avant et aprés décrochage donne la fréquence de décrochage ou plutét, 
comme pour le bruit, le produit de cette fréquence par le nombre de secteurs, 
qui ne peut étre déterminé par cette méthode. 

La figure 12 donne une telle analyse effectuée jusqu’a 500 Hz. La fréquence 
de décrochage est de 14,2 t./s (NV = 1750 t./m soit 29,16 t./s). 

La méthode est encore plus précise que l’analyse du bruit, car le sommet de 
vibrations est plus accentué (figure 10). 

La figure 10 donne le détail de l’analyse au voisinage de la fréquence du 
décrochage, effectué avant et aprés décrochage pour une ouverture déterminée 
et pour quatre vitesses de rotation. 


3.3. Comparaison des diverses mesures 


Les diverses méthodes de mesures que nous avons utilisées (tachymétre a 
lame vibrante, enregistrement de la pression instantanée, analyse spectrale du 
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-e-e-- Analyse du bruit 


—o—o— Analyse des vibrations 
—ss— N=2250% N, = 18,25 ¢,/s 
—e—e— N=2000t, N,=1635c¢/s 
ee 0 N=1750n N=14.3 c/s 
ee NN =1500, p= 123 c/s 
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Aprés decrochage 


Avant décrochage 


8 10 12 14 16 18 [Hz] 
Figure 10 


Analyse harmonique détaillée du bruit et des vibrations avant et aprés décrochage, au voisinage de la 
fréquence du décrochage tournant, pour diverses vitesses de rotation (analyse 4 ouverture constante). 


Figure 11 


Placement du capteur de vibrations sur l’enveloppe de la machine. 
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bruit et des vibrations) donnent des résultats dont la dispersion est extré- 
mement faible. | 


Les valeurs déduites des enregistrements de la pression instantanée doivent 
étre considérées comme exactes. 


Nous avons signalé la précision des analyses spectrales, de l’ordre de 05%. 
La mesure par tachymétre 4 lame vibrante est moins précise, 4 cause de la 
caractéristique de résonance de l'appareil assez aplatie. 


4. Conclusions 


Dans le présent travail, nous avons établi: 
a) Les lois de similitude dans la région de fonctionnement avec décrochage. 
b) La loi de similitude de la vitesse du décrochage tournant brusque. 

Nous avons montré en outre qu'il était possible de déterminer avec une 
bonne précision (-+ 0,5%) la vitesse de rotation du décrochage tournant par 
lanalyse spectrale du bruit ou des vibrations induites par un rotor isolé. 
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Summary 


It is known that when an axial flow compressor running at constant speed is 
throttled down, a phenomenon of rotating stall may occur in one or several inter- 
mediate stages. The velocity of propagation of the latter phenomenon has to be 
investigated as it determines the frequency of the blade vibrations excited by 
this non-steady stall. 

This investigation was carried out on a single rotor. The similarity laws in the 
working region associated with rotating stall were found and the similarity laws 
relative to the velocity of propagation of stall were also determined when a 
discontinuity in the characteristic curve was associated with this phenomenon. 

It is also shown that the velocity of rotating stall may be determined from 
an induced phenomenon: sound or vibration analysis. 


(Regu: le 6 octobre 1956.) 
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Integralgleichungen fiir einige Randwertprobleme 
fir Gebiete mit Ecken 


Von Kurt ARBENZ, Ziirich!) 


Die klassischen Integralgleichungen fiir die konforme Abbildung [1] *) und fir 
Spannungsprobleme haben den Nachteil, dass sie nur fir Gebiete mit stetig 
gekriimmter Randkurve gelten. Es werden neue Integralgleichungen mitgeteilt, 
die auch noch fiir Gebiete mit Ecken gelten und doch numerisch bequem ausge- 
wertet werden kénnen. Zum Beweis der Existenz der Lésungen und zu Aussagen 
iiber die Eigenwerte werden Hilfsmittel der Funktionalanalysis benutzt ilk 


1) Institut fiir angewandte Mathematik der ETH. ae, ; 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 254. 
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1. Konforme Abbildung 


Das einfach zusammenhangende Gebiet G mit stiickweise glattem Rand C 
soll konform auf den Einheitskreis abgebildet werden, so dass der innere Punkt 
Q in den Nullpunkt tibergeht. 

s und ¢ sind die zu den Punkten s,¢ auf C gehdrenden Parameterwerte, 
nachdem auf C ein beliebiger Parameter gewahlt wurde, der nicht die Bogenlange 
zu sein braucht. a ist ein fester Punkt, der auf einem glatten Bogen liegt und 
zugleich Anfangspunkt des Integrationsintervalls ist. 


Figur 1 
Illustration der Verhaltnisse des Gebietes G. 


Fiir den Polarwinkel o(s) des Bildpunktes von s auf C gilt die Integral- 
gleichung 1 
2 w(s) = o(s) — ap) a, Os; t) 2 ola) = 2rea(a): (1) 


Die Integration ist im Sinne von STIELTJES gemeint. Diese Gleichung kann als 

Modifikation der Gerschgorinschen Gleichung angesehen und auch durch Inte- 

gration aus einer von STIEFEL benutzten Gleichung hergeleitet werden. 
Bezeichnet man mit J den Integraloperator 


a= Zh ay Hs, 9 10), 


so kénnen tiber den Operator A = 1 — I folgende Aussagen gemacht werden: 

1. A ist ein linearer und beschrankter Operator im Raume der auf C stetigen 
Funktionen. 

2. Die Eigenwerte 4 von A sind reell und liegen im offenen Intervall 0< 4< 2. 

3. A= 0 ist kein Eigenwert, das heisst, der Operator A ist regular, und somit ist 
die Lésung von (1) eindeutig. A = 2 ist kein Eigenwert. 

4. Die Haufungspunkte der Eigenwerte liegen im abgeschlossenen Intervall 
[1 — da/z, 1+ da/m]. da ist der grésste Sprung, den die Halbtangente in 
Richtung der positiven Orientierung mit einer festen Richtung bildet. 

5. Der Kern der Integraltransformation ist beschrankt, so dass bei der nume- 
rischen Auswertung keine Schwierigkeiten entstehen. 


Die konforme Abbildung im Innern ist bis auf einen multiplikativen Faktor 
durch das Stieltjessche Produkt gegeben 


Zo p 3 
Hey) = 2 TT a ep. 


Das Stieltjessche Produkt ist das Analogon des Stieltjesschen Integrals. 
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2. Das Spannungsproblem einer ebenen Schale 


Die Scheibe G sei Kraften unterworfen, die-in ihrer Ebene liegen und am 
Rand C angreifen. (Uber G gelten die gleichen Voraussetzungen wie in 3) 
t= (t,, t,) sei die als stiickweise stetig vorausgesetzte Belastung auf dem Rand C 
in Funktion der Bogenlange s. 

Bezeichnet man mit U die Airysche Spannungsfunktion, so hat man folgendes 
Randwertproblem zu lésen: 


Eo aa = 0 im Innern: 
s Ss 
a =| t,(o) do = f(s), U, a t,(a) do = g(s) auf C. 
a yi a 


Zur Losung dieser Randwertaufgabe macht man den Ansatz 


UP) = h atx vl P hl) + P aly vA) held (2) 


wobei £, und &, stetige Funktionen auf C sind. Zu ihrer Bestimmung hat man das 
System von Integralgleichungen 


gis) = a hls) + P dy a4(0) halt) + Ff ay 0082 96(t) yl) 


1 . 
ae f d, sin 2 o,(¢) Ro(é) . 


Figur 2 
Darstellung des Winkels ¢,(¢) fiir zwei verschiedene Lagen von ?. 


Diskrete Einzelkrafte 


Die Einzelkraft k? greife im Punkt s, auf C an, dann ist 
s 


HGS = i t,(0)do= — Ski, os) = : 1,(0) do = Sih. 


= 
¢ s~ss a 
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f und g sind Treppenfunktionen. Die Operatoren g und sin[45° + 2 9] sind nur 
fiir stetige Funktionen definiert. Durch das folgende Verfahren lasst sich das 
Problem auf ein solches fiir stetige Randwerte reduzieren. Man macht fiir U 
den Ansatz 


U = U* + SU, [U* wie in (2)] 
a 


mit U; = a; [% 95,(P)] + Oly %5,(P)] - es 
Diese U, dienen zur Hervorbringung der richtigen Singularitaten, wobei die 
Faktoren a; und b,; so bestimmt werden, dass die Funktionen 


stetig auf C sind. 
Nun lést man das Randwertproblem: 


1. JdM4U* = 0 im Innern; 
Be" sav fh LU sare Ah ae 


Bemerkungen 


a) Die Spannungen erhalt man durch Differentiation der Airyschen Spannungs- 
funktion U, wobei die Differentiation unter dem Integral vorgenommen werden 
kann. 

b) Die Kerne im System (3) sind beschrankt. 

c) Um ein Stieltjessches Integral numerisch auszuwerten, muss man auf dem 
Rand eine beliebige Einteilung vornehmen. Die einzige Vorschrift ist, dass jede 
Ecke Teilpunkt ist. 

d) Es wurde ein Dirichletsches Randwertproblem fiir ein Quadrat gelést, fiir 
die numerische Auswertung wurden auf dem Rand 8 Teilpunkte gewahlt. Trotz 
der groben Einteilung iiberstieg der relative Fehler der Funktion U im Innern 
SHY, satel alien ; 


LITERATURVERZEICHNIS 


(1] G. Brrxuorr, D. M. Younc, und E. H. ZARANTONELLO, Effective Conformal 
Transformation of Smooth, Simply Connected Domains. 
[2] F. Riesz und B. Sz.-Nacy, Legons d’analyse fonctionelle. 


Summary 


Concerning conformal mapping and problems of the theory of elasticity, the 
classical integral equations are modified in order to treat the cases of not smooth 


boundary curves. As the kernels are still bounded, no difficulties of numerical 
computation arise. 


(Eingegangen: 22. Dezember 1956.) 
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Varia — Miscellaneous — Divers 


International Congress of Mathematicians, Edinburgh, 
August 14th - 21st, 1958 


At the invitation of the City and University of Edinburgh and the Royal 
Society of London, the International Congress of Mathematicians will meet in 
Edinburgh from August 14th to August 21st, 1958. His Royal Highness the 
Duke of Edinburgh has graciously consented to extend his patronage to the 
Congress. 

The Executive Committee is inviting a number of mathematicians to deliver 
one-hour and half-hour addresses. There will also be daily sessions devoted to 
fifteen-minute communications. Those who wish to present such communications 
will have an opportunity of offering to do so when they receive the Second 
Communication. 

There will be eight sections, namely: 1. Logic and Foundations. 2. Algebra and 
Theory of Numbers. 3. Analysis. 4. Topology. 5. Geometry. 6. Probability and 
Statistics. 7. Applied Mathematics, Mathematical Physics and Numerical Analysis. 
8. History and Education. 

Those who wish to receive further information about the Congress are requested 
to communicate their name and full address as soon as possible to the secretary: 
Dr Frank Smithies, Mathematical Institute, 16 Chambers Street, Edinburgh, 1, 
Scotland. Towards the end of 1957 they will receive a Second Communication 
giving more detailed information and specifying the fees to be paid by members; 
a form of application for membership of the Congress will be included. 

F. Smithies 


Buchbesprechungen — Book Reviews — Notices bibliographiques 


Elementare Differentialgeometrie. Von W. Haack (Birkhauser Verlag, 
Basel und Stuttgart 1955). 239S., 25 Abb.; gebunden Fr. 22.-, broschiert Fr. 18.70. 

Dieses Werk ist aus dem Biichlein Diffeventialgeometrie I desselben Verfassers 
hervorgegangen, das 1947 in der Wolfenbiitteler Verlagsanstalt erschienen ist. 
Sein Inhalt umfasste eine Darstellung der klassischen Kurven- und Flachen- 
theorie auf vektorieller Grundlage. In der nun vorliegenden Elementaren Diffe- 
ventialgeometrie bringt Haack dariiber hinaus eine Einfiihrung in die Cartanschen 
Methoden des «repére mobile» (bewegliches Dreibein) und der schiefen Differen- 
tialformen sowie eine eingehende Behandlung einiger Existenzprobleme der 
Flachentheorie. 

Die Darstellung ist leichtverstandlich gehalten. Haacks Buch erleichtert auf 
diese Weise dem Anfanger den Zugang zu den wichtigsten Ergebnissen und Ar- 
beitsmethoden der Differentialgeometrie. Insbesondere sei auf die Auseinander- 
setzung der Ideen von CarTAN hingewiesen, die jetzt nach und nach in der 
deutschsprachigen Literatur zur Infinitesimalgeometrie Eingang finden. Der Ver- 
fasser entwickelt zuerst die Grundprobleme mit elementaren Hilfsmitteln und 
dann hernach nochmals von neuem im Sinne von Cartan. In den Schlusskapi- 
teln verwendet Haack dann nur noch den schiefen Differentialkalkiil. So wird 


a. 
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der Leser nach und nach mit den modernen Arbeitsmethoden der Differential- 
geometrie vertraut gemacht. Im Hinblick auf die Anwendungen in der Geodasie 
und in der Kartographie sind einige Fragen tiber Flachenabbildungen etwas 
breiter dargestellt. Einzig das Schlusskapitel (Existenzsétze der Flachentheorie) 
verlasst den Rahmen des Elementaren und hat mehr den Charakter einer Mono- 
graphie iiber das Problem der Realisierbarkeit einer vorgegebenen Metrik und 
das Problem der Flachenverbiegung. M. Jeger 


Creep and Fracture of Metals at High Temperatures. Proceedings of 
a Symposium held at the National Physical Laboratory (Her Majesty’s Stationary 
Office, London 1956). 419 S., 258 Fig.; £ 1.10.0. 

Die 22 Vortrage an dieser Arbeitstagung befassten sich mit Arbeitshypothesen 
und Theorien zur Klarung der komplexen physikalischen Vorgange im atomaren 
Gebiet, die bei der plastischen Verformung, der Kristallerholung, der Rekristalli- 
sation, der Alterung und Ausscheidung auftreten und die zum Kriechen und 
Bruch von Einkristallen und polykristallinen Haufwerken fiihren; dies alles stets 
im Hinblick auf die steigenden Anforderungen der Technik an die Dauerwarm- 
festigkeit metallischer Werkstoffe. 

Manche Theorien blieben in der Diskussion umstritten, andere konnten durch 
‘Versuchsergebnisse im Makrogebiet gut unterstiitzt werden, wieder andere for- 
. derten wertvolle Hinweise fiir die Weiterentwicklung kriechfester Legierungen . 
* zutage. 

Die Fiille der Theorien, Versuchsergebnisse und Literaturangaben wurde 
zweckmassig in vier Hauptabschnitte unterteilt: 1. Deformationsvorgange in 
-einfachen Stoffen, 2. Kriechwiderstand komplexer Stoffe, 3. Physikalische Bruch- 
theorien, 4. Vorgange in der dritten Kriechstufe und beim Bruch. 

Allgemein wurden bisherige Theorien tiber die Fortpflanzung linearer, gewen- 
delter, raumlicher usw. Stor-, Locker- und Leerstellen im Atomgitter weiter aus- 
gebaut und damit ein reiches und interessantes Material zur Weiterentwicklung 
der physikalischen Metallkunde vorgelegt. 

Zu wiinschen ware, dass fiir solche und ahnliche junge wissenschaftliche 
Spezialgebiete die Terminologie und Begriffsdefinitionen in den verschiedenen 
Sprachen mit der raschen Ausweitung der Forschung Schritt halten wiirden. 

E. Bickel 


Raum- und Bauakustik fiir Architekten. Von WILLI FuRRER (Birkhauser 
Verlag, Basel und Stuttgart 1956). 200 S., 160 Abb.; Fr./DM 27.50. 

Infolge der neuen Baumethoden und des wachsenden Larmes werden die 
akustischen Probleme immer schwieriger: ihre Lésung macht die Kenntnisse 
theoretischer Grundlagen und experimenteller Messergebnisse notwendig. Das 
kiirzlich erschienene Buch von Prof. FuRRER gibt auf vorztigliche Weise Anleitung 
dazu, angepasst an die Méglichkeiten und Bediirfnisse der Architekten. 

Der erste Teil des Buches behandelt akustische Grundbegriffe, es werden die 
physikalischen Eigenschaften des Schallfeldes und Schallempfindens sowie mess- 
technische Einzelheiten besprochen. 

Der nachste Abschnitt tiber Rawmakustik behandelt die Ausbreitungsverh4lt- 
nisse nach geometrischer, statistischer und wellentheoretischer Methode. 

Im dritten Teil, Bauakustik, werden die notwendigen Kenntnisse iiber die 
- Schallschluckung (Absorption) und Schalldammung (Isolation) vermittelt. 

Der gesamte Stoff ist durch viele Beispiele sehr klar dargestellt. Dank der 
zahlreichen Tabellen und Kurven dient dieses Lehrbuch auch als Nachschlage- 
werk und wird in Fachkreisen gewiss freudige Aufnahme finden. J. Martony 
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Die Transformatoren 


von Dr. oH. c. Dr. TEcHN. ING. MILAN VIDMAR 


Ordentliches Mitglied der Slowenischen Akademie der Wissenschaften und Kiinste 
Ordentlicher Professor an der Universitat Ljubljana 
Direktor des Instituts fiir Elektrizitatswirtschaft in Ljubljana 
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Dritte, vollstandig umgearbeitete Auflage 
630 Seiten mit 321 Abbildungen im Text, davon 2 Falttafeln 
Ganzleinen Fr. 68.- (DM 68.—), broschiert Fr. 64.- (DM 64.-) 
: 1956 


Inhaltsverzeichnis 


I. Das Wesen des Transformators — II. Die Gestalt des Transformators — III. Das Wesen des Drehstrom- 

transformators — IV. Die Warmeabfuhr aus dem arbeitenden Transformator — V. Der Eisenkern. Der 

Leerlauf, Der Einschaltstromstoss — VI. Die Wicklungen — VII. Randprobleme des Transformatorenbaues 
VIII. Der Entwurf. Beispiele von ausgefiihrten Transformatoren. 


Die Transformatorenprobleme haben eine verhaltnismassig gewaltige Vergangenheit hinter sich, sind 
indessen immer noch nicht vollstandig gelést. Es gibt eine Unmenge von mitwirkenden Faktoren, die den 
Bau eines grossen Transformators beeinflussen, und hinter jedem Faktor steht ein Problem. 

«Die Transformatoren», deren dritte Auflage nun herauskommt, haben durch den Erfolg ihrer beiden 
vorangegangenen Auflagen bewiesen, dass sie die Problematik des Transformatorenbaues richtig ange- 
packt und zufriedenstellend behandelt haben. Seit dem Erscheinen der zweiten Auflage sind allerdings 
mehr als 30 Jahre vergangen, so dass die dritte Auflage den Rahmen der zu behandelnden Probleme 
ganz erheblich erweitern musste. Sie hat vor allem die Aufgaben zu behandeln gehabt, die mit den sehr 
hohen Spannungen im Transformatorenbau aufgetreten sind. Dariiber hinaus gab es neue konstruktive 
Probleme, die zum Teil mit den ganz grossen Leistungen zusammenhangen, zum Teil mit dem Auftauchen 
neuer Baustoffe verbunden sind. 

Dies alles soll die dritte Auflage «Der Transformatoren» dem Sonderfack, in das sie ausgerichtet ist, 
bringen. Sie kommt trotz der erweiterten Problematik mit einem etwas bescheideneren Umfang aus als die 
vorangegangene zweite Auflage. Trotzdem bemiiht sie sich, die praktische Anwendung der beschriebenen 
Problemlésungen anhand von Beispielen zu erleichtern. ? 
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